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Preface 


This  document  comprises  a  collection  of  technical  notes  'n  various  aspects  of  high- 
power  microwave  (HPM)  generation,  coupling,  and  interaction.  The  notes  were 
written  between  October  1983  and  February  1986,  and  were  originally  distributed 
rather  informally  as  Microwave  Notes  to  members  of  the  HPM  community.  It  has 
since  seemed  worthwhile  to  edit  and  transcribe  them,  and  to  publish  them  as  a 
single  report.  Time  and  resources  did  not  permit  either  as  detailed  an  editing 
job  as  would  have  been  desirable,  or  recalculation  and  plotting  of  the  numerical 
results  in  a  more  elegant  form.  It  is  to  be  hoped,  however,  that  this  integrated 
document  will  prove  useful  to  the  community. 

The  selection  of  topics  addressed  in  the  Microwave  Notes  was  influenced  by  ques¬ 
tions  which  arose  at  various  stages  of  the  author’s  involvement  in  HPM-related 
activities.  The  subject  of  pulse  radiation  from  aperture  antennas  subject  to  turn¬ 
on  time  errors  ^discussed  ITT  for  example,  arose  in  the  context  of  a  pro¬ 

posed  HPM  generation  scheme.  Consideration  of  the  scattering  cross-section  of  a 
simple  dipole  and  of  its  dual  problem,  the  effective  area  of  a  thin  slot  fdiscnsscd  c 
was  motivated  by  the  need  to  understand  the  penetration  of  such  slot 
apertures  by  HPM  pulses.  The  slot-aperture  penetration  problem  was  considered 
further  in  Note  5,  where  a  monTcompficatecTsTot-to-internal-wire  coupling  config- 
uration  was  addressed,  and  in  Note  7,  where  the  problem  of  coupling  through  a 
deep  slot  was  considered.  The  “front-door”  coupling  problem  was  also  of  interest; 
hence  the  treatments  of  pulse  propagation  in  waveguides  ( Note  3)  and  of  coupling 
to  a  filamentary  load  in  a  waveguide  (Note  6).  Interest  in  estimating  the  HPM 
power  that  might  be  scattered  into  an  antenna  located  over  a  rough  surface  by 
an  incident  field  at  which  the  antenna  was  not,  in  fact,  “looking”  motivated  the 
study  of  that  problem  ( Note  4). 

The  order  of  the  chapters  in  this  report  follows  that  of  the  original  Notes.  It  is 
perhaps  worthwhile  to  point  out  that  the  analyses  presented  herein  by  no  means 
represent  original  research.  They  were  intended  simply  to  make  use  of  tried-and-  j 
true  analytical  methods  for  the  study  of  problems  of  current  interest.  Thus  the 
reader  may  be  disappointed  by  the  lack  of  extensive  references.  Suffice  it  to  say  *4 
that  Harrington’s  Time-Harmonic  Electromagnetic  Fields  and  the  Handbook  of  loa 
Mathematical  Functions  (M.  Abramowitz  and  I.  A.  Stegun,  eds.)  were  most  useful 
during  the  preparation  of  the  Notes.  . — ~ 


«»/ _ 
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It  is  a  pleasure  to  acknowledge  the  contributions  of  Dr.  Eric  P.  Wenaas  of  JAY- 
COR,  who  originally  suggested  that  I  write  the  Microwave  Notes-,  Dr.  Kelvin  S. 
H.  Lee  of  Kaman  Sciences  Corporation,  with  whom  I  have  had  the  pleasure  of 
discussing  the  problems  considered  herein  (as  well  as  many  others);  Mr.  Gary 
Kwitkoski  of  DNA,  who  encouraged  the  preparation  of  this  report;  and  Ms.  Susan 
G.  Brown  of  JAYCOR,  who  transcribed  the  original  Notes  into  TgX  format  and  in¬ 
tegrated  them  into  this  document.  Errors  and  omissions  (which  are  surely  present) 
are  my  responsibility  alone. 
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SECTION  1 

PULSE  RADIATION  FROM  APERTURES 
WITH  RANDOM-PHASE  EXCITATION 


Summary:  We  examine  the  field  radiated  by  a  circular  aperture  antenna  which 
is  excited  by  a  pulsed  field  subject  to  turn-on  time  errors  randomly  distributed 
across  the  aperture.  In  particular,  we  calculate  the  expected  value  and  the  auto¬ 
correlation  of  the  impulse  response  of  the  radiated  field  when  the  turn-on  time  is  a 
random  function  of  location  in  the  aperture.  It  is  assumed  that  this  function  is  a 
stationary  Gaussian  random  process  with  zero  mean  and  an  isotropic  correlation. 
Analytical  and  numerical  results  are  presented  showing  the  pulse  stretching  that 
occurs  as  a  result  of  the  non-zero  antenna  fill  time  for  off-axis  observers  and  of  the 
random  turn-on  time  errors.  We  also  present  results  showing  the  energy  density 
radiated  when  the  aperture  field  is  a  short  cw  pulse  in  the  time  domain. 

In  order  to  assess  the  potential  threat  posed  by  microwave  directed-energy  weapons 
one  must  first  understand  the  characteristics  of  the  signals  which  can  be  radiated 
by  such  weapons.  These  signals  depend  upon  the  microwave  energy  source  and 
upon  the  antenna  system. 

In  this  section  we  examine  the  effect  of  the  antenna  itself  on  the  radiation  of  a 
pulsed  signal.  In  particular,  we  consider  a  circular  aperture  antenna  excited  by  a 
field  which  is  spatially  uniform  except  for  random  variations  caused  by  localized 
“turn-on  time”  errors.  Using  this  model  we  investigate  the  effects  of  the  finite 
antenna  fill  time,  and  of  the  turn-on  time  errors,  on  the  impulse  responses  observed 
in  the  far  field.  Because  of  the  assumed  random  character  of  the  aperture  field,  we 
consider  certain  statistical  descriptors  of  the  far-zone  radiated  field:  the  expected 
value  and  the  autocorrelation  of  the  impulse  response.  These  quantities  can  be 
combined  by  convolution  with  a  given  aperture  field  to  yield  the  expected  signal 
radiated  by  an  aperture  field  of  any  specified  time  dependence. 

Random-phase  effects  have  been  considered  in  the  past  for  cw  radiation  from 
aperture  antennas  (lj.  Such  errors  produce  a  reduction  in  the  on-axis  antenna 
gain  and  a  broadening  of  the  antenna’s  radiation  pattern  lobes.  These  phenomena 
can  be  understood  on  the  basis  of  power  transfer  from  the  “coherent”  or  error-free 
radiation  pattern  into  the  “incoherent”  pattern.  The  effect  depends  in  detail  on 
the  variance  and  the  correlation  length  of  the  phase  errors. 


1 


We  shall  examine  these  effects  in  the  time  domain.  The  general  formulation  of  the 
problem  is  given  next,  followed  by  analytical  and  numerical  results  for  the  field 
and  energy  density  impulse  responses.  Implications  of  these  results  for  microwave 
pulse  radiation  are  discussed. 

1.1  FORMULATION. 

We  consider  a  circular  aperture  antenna  of  radius  a  located  in  the  plane  2  =  0,  as 
shown  in  Figure  1.  A  scalar  field  Exa(p,(f>;t)  is  imposed  over  the  aperture.  This 
aperture  field  is  taken  to  be  of  the  form 

Exa(p,<f>]t)  =  Exa[p,<t>\t  -  T(p,<t>)\  (1) 

where  r(p,  cj>)  is  assumed  to  be  a  two-dimensional  Gaussian  random  process  having 
zero  mean,  variance  rj,  and  isotropic  correlation  function  r(£).  Thus 

<  T(p)  >~  0  (2) 

<  *{pW)  >=  ror(l  P~P'  I)  (3) 

with  the  brackets  denoting  the  expected  value;  p  and  p '  denote  positions  in  the 
aperture. 

In  what  follows  we  shall  assume  that  the  aperture  field  is  spatially  uniform  except 
for  the  variation  in  turn-on  time  r;  thus 


EXa{p,  4>\  t)  =  Eo[t  -  t(p,  <i>))  (4) 

The  Fourier  transform  of  the  aperture  field  is 

Era(P,<f>)  =  I"  E0[t-T(p,<f>)]e-^dt 

J  —CO 

too 

=  /  E0(t)e-]“ldt 

J—oo 

=  (5) 

where  Fc  denotes  the  Fourier  transform  of  the  aperture  field  in  the  absence  of 
turn-on  time  errors.  The  far-zone  radiated  field  corresponding  to  the  assumed 
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aperture  field  is  well  known;  it  is  given  in  the  frequency  domain  by 


Ex(r,0,<j>) 


jkEo(ju>)  ^-jkr  f 2  fa  ejkp' sin  6  cos(<t>-4>')  _ 

2irr  Jo  Jo 


(6) 


g-jWp'^y 

wherein  k  is  the  free-space  propagation  constant1  and  Ex  is  expressed  in  terms  of 
the  spherical  coordinate  system  shown  in  Figure  1. 


Figure  1.  Aperture  geometry. 


In  what  follows  we  shall  make  use  of  the  mean  and  correlation  functions 

<  >=  e_wa^/2  (7) 


(8) 


lk  =  UyJ /jq€q,  where  /i0  and  c0  are  respectively  the  permeability  and  the  permittivity  of  free 
space. 
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wherein  the  function  /(•)  is  defined  by  Equation  (8).  We  shall  calculate  the 
expected  value  and  the  autocorrelation  of  the  radiated  field  under  the  assumption 
that  the  aperture  field  in  the  absence  of  turn-on  time  errors  is  an  impulse  in  the 
time  domain. 


1.2  IMPULSE  RESPONSE  OF  THE  RADIATED  FIELD. 

The  time-domain  radiated  field  is  given  by 

Ex(r,6,<t>;t)  =  —  j  Ex(r,  0,  <f>)ejutdu} 

Z7T  J—o o 

and  the  expected  value  of  this  field  is 

1 

<  Ex(r,  0,  fat)>=  —  <  Ex(r,  0 ,  4>)  >  e]UJtdu> 

Z7T  J  —  co 


(9) 


(10) 


by  virtue  of  the  linearity  of  the  integration  and  expectation  operations.  Using 
Equations  (6)  and  (7),  we  obtain 


<  Ex(r,0,<t>)>  =  /2e-jkr 


2i rr 

2  7T  ra 


(11) 


na 

. 


,jkp'  s\n0  cos^-t') p, dp'  dtf 


The  integral  in  Equation  (11)  can  be  readily  expressed  in  terms  of  a  Bessel  function 
of  the  first  kind,  yielding 


< 


Ex(r,0,<j>)  >=  j-afe(JU?) e~Jtro/2e~jkTJl(ka sin 9) 
rsmv 


We  shall  find  it  useful  to  employ  the  relation 

■/,(*)  =  -  /' 

J7T  J- 1 


1  i««  Ud%i 


\/l  —  U 2 


so  that 


<  Ex{r,0,<f>)  >=  ^Me-}kre-u,lTj/2  p  ejkaus 

7rrsin0  J- 1  y/l  —  u 2 


(12) 


(13) 


(14) 
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Now  if  E0(t)  is  a  delta  function,  Eo{jw)  =  1  and  we  define 


<  h(t  -  r/c)  >  =  — [°°  eMt~r/c)e~^To/2du;  ■  (15) 

2w2r  sin  u  J—oo 


cjka u  sing  U^U 


as  the  expected  value  of  the  impulse  response  of  the  radiated  field.  Suppressing 
the  time  delay  r/c,  interchanging  the  order  of  integration,  and  performing  the 
integration  over  w,  we  obtain 


<  h(t)  >= 


ntjr  sin# 


£  4v;; 


udu 


\J  1  —  u2 


(16) 


wherein  we  have  defined  the  Gaussian  function 

G(x\ -  ~%=~e~*i,2a*  C7) 

V  27T (7 

and 

t f  =  -  sin  9  (18) 

c 

denotes  the  radial  “fill  time”  of  the  aperture  for  an  observer  at  an  angle  0  with 
respect  to  the  z-axis;  c  =  1  / y/no^o  denotes  the  speed  of  light  in  free  space. 

Now  it  is  easy  to  show  that  h(t)  reduces,  as  it  should,  in  the  limit  tq  — ♦  0  to 


lim  <  h(t )  >=  — - — : — -  - - - 

ro-°  irrt fsm0  (</</)2 


(19) 


and  also  that  <  h(t)  >  is  an  odd  function  of  t  for  any  value  of  the  parameter 
r0/tf.2  We  have  calculated  the  normalized  impulse  response 


m  = 


nrtj  sin  0 
a 


<  h(t)  > 


(20) 


as  a  function  of  normalized  time  t/tj,  for  various  values  of  the  parameter  T0/tj. 
The  results  are  shown  in  Figure  2.  It  is  evident  from  the  curves  shown  that 

2We  remark  that  <  h(t)  >  is  in  general  not  causal.  This  is  a  consequence  of  our  assumption 
that  the  turn-on  time  is  a  Gaussian  random  process. 
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the  effect  of  increasing  the  value  of  r0  (the  standard  deviation  of  the  turn-on 
time)  smooths  and  broadens  the  expected  value  of  the  impulse  response.  It  is  not 
difficult  to  show  that  when  To/tf  becomes  comparable  to  or  exceeds  unity, 


so  that  the 
than  by  tf. 
for  the  case 


— irrtjsmO 
a 


<  h(t)  >=  - 


(21) 


width  of  the  impulse  response  is  ultimately  controlled  by  r0,  rather 
We  show  in  Figure  3  a  plot  of  ( To/tf)2f{t )  vs.  normalized  time  f/r0 
where  r0  dominates  the  impulse  response. 


Figure  2.  Normalized  impulse  response  f{t)  vs.  normalized  time  t/tf,  for  r0/tj= 
0.1,  0.5,  and  1.0,  0  <  t/tf  <  =  -/(<)• 


It  is  evident  that  the  effects  of  the  random  turn-on  time  are  dominant  for  small 
observation  angles  9.  If  we  take  as  a  maximum  of  tf  the  value  a/c ,  we  see  that 
turn-on  time  effects  will  dominate  the  impulse  response  for  all  observation  angles 
if  cr0/a  is  comparable  to  or  greater  than  unity. 

Since  the  impulse  response,  whether  it  is  dominated  by  tf  or  by  r0,  has  a  double 
peak,  the  radiated  field  will  exhibit  the  same  feature  if  the  aperture  field  has  a 
duration  short  compared  to  that  of  the  impulse  response.  Thus  a  single  short 


6 


pulse  in  the  aperture  can  appear  in  the  radiated  field  as  a  double  pulse  under 
certain  conditions.  The  two  components  of  the  double  pulse  will  be  separated 
by  a  time  of  order  2 tf  or  2 r0,  whichever  is  larger,  and  they  will  be  opposite  in 
sign.  Under  these  conditions  the  antenna  cannot  be  considered  a  “faithful”  pulse 
radiator,  because  of  the  breakup  of  the  radiated  time-domain  signal. 

The  expected  value  of  the  impulse  response  does  not  involve  the  spatial  correlation 
of  the  turn-on  time  over  the  aperture.  To  show  the  effect  of  this  descriptor  of  the 
aperture  fields,  we  next  consider  the  autocorrelation  of  the  radiated  field  impulse 
response. 


Figure  3.  Normalized  impulse  response  (toA/)2 /(0  vs  normalized  time  t/r0  for 
the  case  where  r0  dominates  the  response;  f(—t )  =  —/(<). 
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1.3  AUTOCORRELATION  OF  THE  RADIATED  FIELD 
IMPULSE  RESPONSE. 


The  expected  value  of  the  energy  spectral  density  in  the  radiated  field  is 

U(u)  =  ±-<\E,(r,$,*)  |2> 

6 0 


(22) 


where  Z0  =  \j  p o/co  denotes  the  intrinsic  impedance  of  free  space.  Using  equations 
(6)  and  (8),  we  have 


U{u) 


k2  1  Eojju)  \2 
4t r2Z0r2 


ra  r2ir  fa  f2it 

/  p'dp'  /  d<f>'  /  p'dp"  f  dcj>" 

Jo  Jo  Jo  Jo 


ejk  sin  6[p'  cos (4>-<t>')-p"  cos . 


(23) 


+  (l  -e~^)f{\p'-p" 

Now  U(u>)  as  given  by  Equation  (23)  above  can  be  expressed  as 

lW|2«2.-«  ,2, 


U{u)  = 


-e  U^T°  ./12(fcasin0) 


Z0r2  sin2  0 


(24) 


where  f  denotes  the  unit  vector  in  the  radial  direction.  If,  as  we  shall  assume,  the 
correlation  distance  (i.e.,  the  distance  over  which  /(|  p '  —  p"  |)  is  close  to  unity) 
is  small  in  comparison  to  the  aperture  radius  a,  the  integral  in  Equation  (24)  may 


8 


be  approximately  evaluated,  yielding  for  U(ui )3 


Zor''  sin  0 


kW  |  I2 


(l  —  e  “2to)  J  f  [u)Jo{ku  sin  0)udu 


The  autocorrelation  of  the  radiated  field  impulse  response,  //(£),  is  defined  as  the 
inverse  Fourier  transform  of  ZqU(u>)/  j  Eo(ju)  |2: 


H{t)  = 


r2  sin2  6 


T~x  {e'^oJ2  {ka  sin  9) 


-f  ^(&/sin#)2(l  —  e~w*T°)  J  } (l v)J0(klv  sin  0)vdv^ 

where  l  denotes  the  characteristic  correlation  length.  H{t)  thus  comprises  the 
sum  of  the  autocorrelation  of  the  expected  value  of  the  impulse  response  and  an 
additional  term  which  depends  upon  the  spatial  correlation  of  the  aperture  field 
turn-on  time.  To  illustrate  the  behavior  of  the  added  term,  consider  the  simple 
case  where  the  aperture  field  turn-on  times  are  perfectly  correlated  (r  =  1)  over 
circular  regions  of  radius  /  (/  <<  a)  and  uncorrelated  outside  these  regions.  Then 
defining 

Fx{t)  =  {e-^JftJfcasinfl)}  (27) 


W)  =  ^-‘{^(Wsin^l-e-^o). 

J  f(lv)J0(klv  sin  0)vdv^ 


3In  the  event  that  the  phase  correlation  distance  is  comparable  to  or  larger  than  the  aperture 
radius,  U(u)  reduces  to  the  simple  result 

r'r/  \  I  To {jw)  i 2  a2  ,2  ,,  •  m 

T(w)  =  J\  (kasmO) 

Zor^sm  9 

which  is  not  dependent  on  the  statistical  properties  of  the  turn-on  time. 
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so  that 


we  find  4 


Fi(t)  =<  A(f)  >  *  <  h(t)  >  (30) 

F2(t)  =  T~x  | ^A:/ sin 0(1  -  e-^rl)Jx{kl sin0)|  (31) 

When  the  correlation  distance  l  is  small  (more  precisely,  if  A;/sin0  <<  1  over  the 
frequency  range  of  interest), 


F2(t)  =  T 


-  {(l*isine)2 (* _e"“!r°!)} 

=  _  (h”£y  £.  [w _ 

\  2  c  )  dt 2  [  2y/n  t0 


Furthermore,  if  r0  is  comparable  to  or  larger  than  tj, 

FJt)  a  -  ( if™*)’  ( — 1  e-’/W 

\  2c  )  dt 2  \  2y/ir  T0 


so  that 


m  3  ~{0£  [(>-£)  (•)*>]  <34> 


Thus  the  fact  that  the  aperture  field  turn-on  times  are  not  highly  correlated  does 
not  lengthen  the  duration  of  the  autocorrelation  function  H(t),  and  we  conclude 
that 

1.  If  the  aperture  field  turn-on  times  are  highly  correlated,  the  autocorrelation 
duration  is  controlled  by  the  fill  time  tj,  essentially  independent  of  tq. 


4The  symbol  *  denotes  the  correlation  operation. 
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2.  If  the  aperture  field  turn-on  time  is  not  highly  correlated,  the  autocorrela¬ 
tion  duration  is  controlled  by  the  larger  of  the  fill  time  tj  or  the  standard 
deviation  Tq. 


In  the  next  section,  we  shall  consider  the  effects  of  the  fill  time  and  the  turn-on 
time  errors  on  the  radiated  field  when  the  aperture  field  is  a  cw  pulse  and  the 
correlation  function  of  the  turn-on  time  is  an  exponential  function: 

r(u)  =  exp  (— iz2//2)  (35) 


The  function  U(u)  for  this  case  is 
I  E0{ju)  |2  a2 


U{u)  = 


+  e 


Z0r 2  sin2  0 
-J2 


{e_a^ro  J2  (fca  sin  0) 


t2  ( kl  .  fl\\cM“ro)2” 
o  I  —  sin  0  )  V  -  7  -7-  exp 

\  2  )  tx  n{n\) 


(36) 


1.4  PULSED-CW  RADIATED  ENERGY  DENSITY. 


We  have  shown  in  the  preceding  sections  that  the  impulse  response  duration  is 
essentially  the  sum  of  the  aperture  fill  time  tj  =  (a/c)  sin  0  and  the  standard 
deviation  of  the  turn-on  time  error  r0.  If  the  aperture  field  is  a  pulsed  signal 
of  duration  2 T,  then  the  radiated  field  will  be  a  “stretched”  pulse  of  duration 
2T  +  2tj  -f  2 r0.  In  this  section  we  consider  the  radiated  energy  density  (J/m2) 
when  the  aperture  field  is  a  Gaussian  pulse  of  the  form 

E0(t)  =  Epe~t2/2t°  cosu >0t  (37) 


Ep  denotes  the  peak  electric  field  strength  in  the  aperture  and  t0  is  a  measure  of 
pulse  width.  The  center  frequency  of  the  pulse  is  u>0- 

The  total  radiated-field  energy  density  is  readily  shown  to  be 


U  = 


Tr^pfoQ2  f°° 

Z0r 2  sin2  0  J- 00 


e-(w- 


[Fi(u>)  +  F2(u;)]  du 


(38) 


with 


Am 


Fi(lo)  =  e  J^(kasmO) 


(u  .  . 

Vf 

Mo)2” 

1  | 

(ki .  y 

(? 8ln  e 

/  n=l 

(  eXP 

n(n!) 

n 

(39) 

(40) 


We  shall  consider  the  most  important  special  case,  that  in  which  to  >>  tj  and 
t0  »  r0,  so  that  the  radiated  pulse  is  a  quasi-cw  signal.  The  total  energy  radiated 
is  readily  found  to  be 


U 


U0 


,2-2 

0T0 


Jl(uQtj) 
M  t/)2 


(41) 


+ 


oo 


£ 


MM2" 

n(n!) 


exp 


where 

n^kyElto 

°  4Z„r2 

denotes  the  energy  density  radiated  in  the  direction  0  =  0  when  there  is  no  turn-on 
time  error. 


Curves  of  U/U0  as  a  function  of  u0tf  are  given  in  Figures  4-22  for  various  values 
of  u0 To  and  l/a.  We  observe  that  as  l/a  (the  normalized  coherence  length  in  the 
aperture)  decreases,  the  effect  of  increasing  the  normalized  standard  deviation 
uj0t0  becomes  more  pronounced,  appearing  first  at  smaller  values  of  (recall 
that  turn-on  time  errors  will  dominate  the  impulse  response  if  r0  >  tj). 
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Figure  8. 


Figure  9. 
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Figure  22.  Normalized  energy  density  U jU o  vs.  :  l/a  —  0.1,  ujqt0  =  7r. 

1.5  CONCLUDING  REMARKS. 

We  have  considered  the  electromagnetic  field  radiated  from  a  circular  aperture 
excited  by  a  field  which  is  uniform  except  for  randomly  distributed  turn-on  time 
errors.  We  have  found  the  expected  value  of  the  impulse  response  and  of  the 
autocorrelation  of  the  radiated  field,  and  have  shown  that  the  shape  and  the 
duration  of  these  functions  depend  upon  the  aperture  fill  time  tf  =  (a/c)  sin  0 
and  the  standard  deviation  of  the  turn-on  time  r0.  The  fact  that  these  functions 
possess  two  peaks  separated  by  twice  the  sum  of  r0  and  t /  means  that  the  field 
radiated  by  a  sufficiently  short  pulsed  aperture  field  will  comprise  two  more  or  less 
distinct  components  separated  in  time.  Thus  a  pulse  radiated  from  the  aperture 
will  “hang  together”  only  if  it  is  at  least  comparable  in  length  to  if  and/or  r0. 

We  have  presented  extensive  numerical  data  for  the  pulse  energy  density  radiated 
by  a  quasi-cw  aperture-field  pulse  and  have  shown  the  influence  of  the  standard 
deviation  and  of  the  correlation  function  of  the  turn-on  time.  In  particular,  for 
a  fixed  value  of  r0,  decreasing  the  correlation  distance  causes  a  reduction  in  the 
energy  radiated  in  the  axial  direction  and  an  increase  in  that  radiated  in  other 
directions:  the  source  becomes,  in  effect,  more  nearly  isotropic.  A  simular  effect 
"i  seen  as  the  standard  deviation  r0  is  increased  for  a  fixed  correlation  distance. 

We  also  note  that  the  effect  of  the  nonuniformity  in  turn-on  time  can  be  quite 
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severe:  the  energy  radiated  in  the  axial  direction  decreases  as  exp(—  UqTq),  where 
u>0  is  the  pulse  center  frequency.  Thus  if  u>0r0  =  7r  (corresponding  to  one-half  cycle 
at  frequency  uq)  the  axial  energy  is  reduced  by  over  40  dB  if  the  turn-on  time 
errors  are  not  highly  correlated  over  the  aperture. 


The  value  of  uj0t  f  associated  with  the  edge  of  the  main  lobe  of  a  circular  aperture 
antenna  pattern  in  the  absence  of  turn-on  time  errors  is  3.832,  the  first  zero  of 
Ji(-).  For  a  short  cw  pulse  to  be  radiated  essentially  as  a  single  pulse  over  the 
corresponding  angular  range  requires  that  ojoto  »  3.832.  If,  for  example,  one 
chooses  co0t0  equal  to  I2n  to  satisfy  this  condition,  then  the  aperture-field  pulse 
contains  six  cycles.  Thus  an  aperture-field  pulse  of  few  cycles  will  be  radiated  as 
a  single  pulse  over  the  (cw)  main  antenna-pattern  lobe.  For  angles  greater  than 
0m  =  sin-1  (ct0/a)  the  radiated  pulse  will  be  broken  into  two  distinct  pulses,  5 
each  of  which  will  have  a  much  lower  amplitude  than  the  single  pulse  radiated  in 
the  main  lobe.  The  transition  occurs  over  the  angular  range 


sin 


-l 


'3.832  /dflV 
.  Vo/. 


<  9  <  sin"1 


(43) 


The  implications  of  these  results  for  microwave  weapon  threat  analysis  are  that 

•  the  principal  threat  comes  from  the  main  lobe  of  the  transmitting  antenna 
pattern  and  comprises  a  quasi-cw  pulse  which  is  stretched  by  an  amount 
depending  on  the  antenna  fill  time  and  the  standard  deviation  of  the  turn¬ 
on  time  error.  The  energy  density  carried  by  this  pulse  can  be  estimated 
from  results  provided  herein. 

•  the  signal  radiated  well  outside  the  main  lobe  of  the  transmitting  antenna 
can  be  expected  to  be  a  “double  pulse”  of  low  amplitude. 

•  turn-on  time  errors  in  the  aperture  field  of  the  transmitting  antenna  can 
drastically  reduce  the  radiated  energy  density. 


5Note  that  if  to  >  a/c,6m  is  a  complex  angle  whose  real  part  is  ff/2;  thus  a  single  pulse  will 
be  observed  at  any  real  angle  9  if  this  condition  is  met. 
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SECTION  2 

SCATTERING  CROSS-SECTION  OF  A  DIPOLE: 
IMPLICATIONS  FOR  MICROWAVE  COUPLING 


Summary:  We  calculate  the  power  scattering  cross-section  of  an  unloaded  dipole 
antenna  and  estimate  the  equivalent  source  strength  of  a  long  thin  slot  ‘Huminated 
by  a  plane-wave  field.  We  also  determine  the  impulse  response  of  the  s  n. tiered 
field  in  order  to  find  the  amount  of  pulse  stretching  caused  by  the  scattering 
process.  It  is  shown  that  under  pulsed  cw  illumination,  the  peak  power  scattered 
by  the  wire  or  transmitted  through  the  slot  can  be  estimated  on  the  basis  of  the 
cw  cross-section  in  cases  of  practical  interest. 

A  fundamental  problem  in  the  study  of  microwave  interactions  with  systems  is 
that  of  determining  the  amount  of  energy  which  could  be  coupled  from  the  exterior 
to  the  interior  of  a  shielded  system  by  means  of  penetration  through  apertures  in 
the  “skin”.  In  this  section  we  investigate  coupling  through  a  narrow  slot  aperture, 
which  may  be  of  resonant  length  or  longer,  in  a  conducting  surface.  Our  approach 
is  an  indirect  one,  in  which  we  first  treat  the  dual  problem  of  scattering  by  a 
straight  wire  in  free  space  and  determine  its  equivalent  power  scattering  cross- 
section  for  continuous-wave  (cw)  illumination.  It  is  then  a  simple  matter  to  obtain 
the  equivalent  transmission  area  of  the  thin  slot. 

The  problem  of  electromagnetic- wave  scattering  by  a  wire  has  received  a  great  deal 
of  attention  in  the  literature.  We  have  elected  to  employ  a  primarily  analytical, 
rather  than  numerical,  approach  to  the  problem  based  upon  the  Hallen  integral 
equation  for  the  wire  current  [2,3].  This  approach  has  been  used  by  Marin  [4] 
to  study  the  natural  modes  of  various  thin-wire  structures  and  by  Bedrosian  [5] 
to  find  the  natural  resonances  and  modes  of  stick-model  aircraft.  In  the  present 
case,  it  leads  to  an  (almost)  closed-form  expression  for  the  scattering  cross-section 
convenient  for  analytical  and  numerical  exploration.  This  expression  then  can  be 
carried  over  with  minor  modification  to  describe  the  electromagnetic  penetration 
through  a  thin  slot. 

Our  principal  interest  is  actually  in  the  pulsed  cw  response  of  the  wire  or  slot. 
In  order  to  estimate  the  total  pulse  energy  scattered,  we  determine  the  impulse 
response  of  the  scatterer.  The  temporal  duration  of  this  impulse  response  defines 
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the  amount  of  pulse  stretching  which  occurs  in  the  scattering  process.  The  peak 
power  scattered  can  be  estimated  from  the  cw  response  of  the  scatterer  when  the 
illuminating  pulse  contains  at  least  a  few  cycles  of  oscillation. 


2.1  DIPOLE  SCATTERING  CROSS-SECTION. 


The  geometry  of  the  dipole  scattering  problem  is  shown  in  Figure  23.  A  thin 
cylindrical  conductor  of  radius  a  lies  along  the  2-axis  between  z  =  ±1/2.  The 
conductor’s  radius  is  assumed  to  be  small  in  comparison  to  its  length  /  and  to  the 
free-space  wavelength  A.  The  scatterer  is  illuminated  by  a  plane  electromagnetic 
wave  polarized  in  the  0-direction  and  incident  at  an  angle  9  with  respect  to  the 
2-axis. 

Because  the  conductor  is  thin,  the  induced  current  and  charge  densities  on  it  are 
essentially  independent  of  angular  position;  and  the  current  is  predominantly  in 
the  axial  direction.  Thus  the  electromagnetic  quantities  of  interest  on  the  scatterer 
are  its  net  axial  current  7(2)  and  line  charge  density  pi(z );  these  quantities  are 
related  by  the  continuity  equation 

dl  ■ 

—  ±  jup(  =  0 
dz 

We  assume  the  time  dependence  exp  (jut)  for  all  field  quantities. 

The  total  power  scattered  by  the  cylindrical  conductor  can  be  expressed  as  [6] 


(44) 


Ps  =  W  f  f  ~r(z)I(z)  -  ^-pi  ‘(z)pi(z') 
-  — 1/2  •'-1/2  " 


sin  k  ]  2  —  z' 
47T  |  2  —  2/  I 


dzdz'  (45) 


where  po  and  c0  denote,  respectively,  the  permeability  and  permittivity  of  free 
space  and  k  =  Uy/poto-  The  incident  power  density  is  expressed  in  terms  of 
the  incident  electric  field  amplitude  Ei  as 

c  -[Ill 

2Z0 


where  Zq  —  \J po/to  denotes  the  instrinsic  impedance  of  free  space.  The  scattering 
cross-section  is  simply 

=  (47) 
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Figure  23. 


Geometry  of  the  dipole  scattering  problem.  The  radius  of  the  dipole 
is  a  and  the  incident  electric  field  is  in  the  0  direction. 


Therefore,  anticipating  that  I(z )  must  take  the  form 

/(*>  =  §/(**) 


(48) 


where  Z'  is  an  impedance  per  unit  length  and  /(•)  is  a  dimensionless  function  of 
kz,  we  find  with  Eqs.  (44)  -  (47)  that 


Ae  = 


Zo 


4  7T  |  Z' 


W/2 

kl/2  ‘ 

''-*1/2  “ 

-M/2 

/•(«)/(*)  - 


df(u)  df(v) 


du  dv 


sin(u  —  v) 
u  —  v 


dudv 


The  problem  of  determining  the  scattering  cross-section  Ae  thus  reduces  to 
of  finding  Z'  and  /(•). 


(49) 

that 


The  current  induced  on  the  scatterer  is  given  approximately  by  the  solution  of 
the  differential  equation  [3] 


^  S infe*”* 

dz2  il 


(50) 
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subject  to  the  end  conditions  7(±//2)  =  0.  The  “antenna  parameter”  Q  is  defined 


ft  =  21n(//a)  (51) 

It  is  an  elementary  exercise  to  show  that  the  solution  of  Eq.  (50)  which  satisfies 
the  end  conditions  is 

II  \  = 

ZoDsin  0{kl  sin  kl) 


kl 


sin  kl  eikzcose  —  2 j  cos  —  sin 

J  2 


•  (kl 

in  I  —  cos  6  j 


sin  kz(  52) 


„  .  kl  (kl  A 

—2  sin  —  cos  —  cos  0  cos  kz 

2  \  2  /  J 


The  function  in  square  brackets,  which  we  shall  identify  as  f(kz),  can  also  be 
written 


/(it)  =  sin  kl  e^ucoa9  —  sin 


kl 

+  cos 0 )  |  eJU 


(53) 


—  sin 


kl, 

yU  -  cos  9) 


It  is  evident  that  I{z)  as  given  in  Eq.  (52)  is  singular  at  the  antenna  resonances 
kl  =  nir  (|  n  |>  1).  To  remove  these  unphysical  singularities,  we  take  the  radiation 
damping  into  account  by  writing  the  factor  (fc/sin  kl)~l  in  the  form  of  a  Mittag- 
Leffler  expansion  [7]  and  then  shifting  its  poles  into  the  upper  half  of  the  complex 
fc-plane.  Thus 

(54) 

where 


1 

D{kl) 


1 

(k‘)2 


1 

+  6  + 


£' 


n=— oo 


(-r 

nir 


1 

kl  —  nir  —  jan 


+ 


1 

nn  +  jan 


(55) 


The  prime  (')  on  the  summation  indicates  that  the  term  for  n  =  0  is  to  be 
omitted.  The  constants  an  are  calculated  using  a  perturbation  analysis  [4]  and 
can  be  written  in  closed  form  as 

an  =  fj[fr  +  ln(2M*)  -Ci(2|n|jr)]  (56) 
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wherein  7  =  0.57721...  is  Euler’s  constant  and  Ci(-)  denotes  the  cosine  integral. 
Using  Eqs.  (48)  and  (52)  with  (54),  we  find  that 

1  —4irjl 


Z '  ZoQsinOD(kl) 

so  the  scattering  cross-section  Ae  is  written 
4 tt/2  t"'‘  t 


(57) 


A‘  fia  sin3  0  |  D(kl)  pLwJ.u„  L 


,  dm 


sin(u  —  u) 
u  —  v 


dudv  (58) 


where  /(«)  is  defined  in  Eq.  (53)  and  D(kl)  in  Eq.  (54).  The  evaluation  of  the 
integral  in  Eq.  (58)  is  tedious  but  straightforward.  The  real  part  of  the  function 
in  square  brackets  can  be  written  as 


n=l 


where 


&[•]  =  a"  cos  (a«(u  ~  V)1 cos  [A*(u  +  u)l 

ai  =  sin2  kl  sin2  0 


(59) 


kl. 

kl. 

a2  =  4  sin 

y(l  +COS0) 

sin 

-(l-cosO) 

<13  =  —2  sin  kl  sin 

a4  =  —2  sin  kl  sin 

a-i  =  cos  0 
a2  —  0 


y(l  +  cos0) 

y(l  -COS0) 

A  =  o 

&  =  i 


(1  —  cos0) 

(1  -(-  COS0) 


(60) 


a3  =  o(l  +  COS  0)  /?3  =  |(1  -  COS0) 

a4  =  |(1  -  cos0)  =  |(1  -fcos0) 


An  integral  of  the  form 
r«/*  -«/* 
-kl/ J  •'-kl/2 


km  f»n  sin(u  —  v)  . 

In  =  /  /  cos  [an(u  -  v)J  cos  [pn(u  +  u)J - auat; 

•'-kin  •'-kin  u  —  v 


(61) 
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can  be  expressed  as 

Jrkl  gjn  (  fkl—£ 

I  — -iCOSQn /  cos  finr)dr)  (62) 

o  £  Jo 

by  means  of  a  straightforward  change  of  variables.  The  latter  form  for  /„  is 
evaluated  as 

/,  =  ^  sin {Si  [(1 +»„-&)*(] 

+Si  [(1  -  an  -  pn)kl\  +  Si  [(1  -  an  +  /3n)kt\ 

+  Si[(l  +  an +  /?„)*/]} 

~^j~cos0nM  {In  [(i  +pn~an)kl]  -Ci[(l  +  J3„  -  an)kl]  (63) 

-  In  [(1  -  an  -  pn)kl\  +  Ci  [(1  -  a„  -  Pn)kl\ 

-  In  [(1  +  an  -  pn)kl\  +  Ci  [(1  +  an  -  0n)kt\ 

+  In  [(1  +a„  +  pn)kl }  -  Ci  [(1  +  an  +  Pn)kl]} 

When  P„  =  0,  we  have 

ln(Pn  =  0)  =  kl  {Si  [(1  +  an)kl\  +  Si  [(1  —  an)fc/]}  (64) 

f  1  —  cos  [(1  •+  an)^  1  —  cos  [( 1  —  an)Ar/]  1 
1  l  +  an  +  1  -an  J 

Thus, 

7i  =  kl  {Si[(l  +  cosfl)fc/]  +  Si[(l  —  cos0)fc/]}  (65) 

f  1  —  cos[(l  +  cos  0)kl\  1  —  cos[(l  —  cos  9)kl\  1 

\  1  +  cos  9  +  1  —  cos  6  J 

I2  =  sin  kl  S\(2kl)  (66) 

-  cos*/(a  +  ln(2*/)-Ci(2*/)] 
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sin  ^(1  —  cos  0)1 

- {Si  [(1  +  cos  9)kl]  +  Si  [(1  -  cos  0)kl]  +  Si(2Jk/)}  (67) 

( 1  —  COS  0) 

COS  (1  —  COS  0)1 

- 7- - — — -  {ln[(l  —  cos0)kl ]  —  Cif(  1  —  cos0)kl] 

(1—  COS0) 

In  [(1  +  cos  0)kl]  4-  Ci[(l  +  cos  0)kl ]  +  ln(2A:/)  —  Ci(2A4)  4-  7} 


sin  f^(l  4-  cos#)] 

I4  =  - ^ - — -^-{S\[(l+coSO)kl]+Si[(l-coS0)kl}  +  Si(2kl)}(68) 

( 1  +  cos  0) 


cos  [y(1  4-  COS0)] 
(1  4-  cos0) 


{In  [(1  +cos0)fc/]  —  Ci[(l  4-cos0)fc/] 


—  In  [(1  —  cos  0)kl\  +  Ci  [(1  —  cos0)fc/]  4-  ln(2fc/)  —  Ci(2 kl)  -(-  7} 
Assembling  these  results,  we  find  that  Ae  is  given  by 

4  id2 


where 


Ae  ft2  |  D{kl)  |2  sin2  0  ^  +  +  ^  + 


f\  =  2  [1  —  cos  kl  cos  (kl  cos  0)]  [7  -f  ln(2fc/)  —  Ci(2fc/)] 


fi  =  2  sin  kl  s'm(kl  cos  0) 


1 2  log  ^tan  0 


4-  Ci  [( 1  4-  cos  0)kl\  —  Ci  [( 1  —  cos  0)kl]  j 

(71) 


/3  =  | A;/ sin2  0sin2  kl  —  2  sin  Ar/[cos(fc/cos0)  —  cos  kl]} 
{Si[(l  +  cos0)A4]  +  Si[(l  —  cos0)fc/]} 


f4  =  —  2  sin2  0  sin2  kl 


sin2  y(l+cos0)]  sin2  y(l  —  cos#)] 

*  A  1  <  /1 


|  sin2  [1 


cos0 


1  —  cos  0 


When  0  =  tt/2,  i.e.,  when  the  cylinder  is  illuminated  at  broadside,  Ae  simplifies 
to 

Airp 

At  =  n2 1  D{ki)  |2^fc/)  (74) 
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with 


g{kl)  =  2(1  —  cos  fc/)(7  +  ln(2fc/)  —  Ci(2&/))  (75) 

+  2  [fci  sin2  kl  —  2  sin  kl(  1  —  cos  fc/)j  Si (kl) 

—  4  sin2  kl  sin2 

2 


We  have  calculated  the  normalized  scattering  cross-section  AJ l2  as  &  function  of 
kl  for  various  values  of  fl  and  9,  and  as  a  function  of  9  with  fi  and  kl  as  parameters. 
The  results  are  shown  in  Figures  24  -  40. 

The  characteristic  resonances  are  apparent  in  the  plots  of  Ae/l 2  vs.  kl,  which 
also  indicate  the  maximum  value  (approximately  2.1)  attained  by  this  quantity  at 
the  fundamental  (half-wave)  resonance.  Thus  the  effective  area  of  the  unloaded 
dipole  near  its  half-wave  resonance  is  of  the  order  of  the  square  of  its  length, 
almost  independent  of  its  radius.  The  effect  of  changing  the  radius,  and  thus  the 
antenna  parameter  ff,  is  to  alter  the  “Q”  of  the  resonance. 


Figure  24.  Ae/l 2  vs.  kl ;  9  =  30°,  fi  =  5. 
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Figure  28.  Ae/l2  vs.  kl ;  9  —  60°,  H  =  10. 
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Figure  29.  AJt2  vs.  ki,  0  =  90°,  ft  =  10. 

The  plots  of  Aji2  vs.  ft  in  Figures  30  -  35  also  show  the  resonant  behavior  of 
the  structure  and  reproduce  the  features  of  the  well-known  radiation  patterns  of 
long-wire  antennas.  In  general,  the  cross-section  tends  to  decrease  as  frequency  is 
increased  from  the  fundamental  resonance.  The  curves  in  Figures  30  -  35  are  for 
cases  where  the  antenna  is  an  integral  number  of  half-wavelengths  long.  When 
the  frequency  lies  between  the  values  shown,  the  curves  are  smoother.  Examples 
are  shown  in  Figures  36  -  40. 
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Figure  32.  Ae/l 2  vs.  6 ;  kl  —  3x,  fl  =  10. 
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Figure  39.  Ae/l2  vs.  0;  kl  =  9ir/2,  Q,  =  10. 
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Figure  40.  Ae/l 2  vs.  0;  kl  =  ll7r/2,  fi  =  10. 


2.2  EQUIVALENT  TRANSMISSION  AREA  OF  A  SLOT. 

The  scattering  cross-section  which  was  derived  in  the  preceding  section  can  be 
used  to  define  an  equivalent  (power)  source  strength  for  a  dipole.  The  equivalent 
transmission  area  of  a  thin  slot  can  be  readily  obtained  from  the  dipole  scattering 
cross-sections. 

The  geometry  of  the  slot-scattering  problem  is  shown  in  Figure  41.  The  slot  has 
width  w  and  extends  from  z  —  —1/2  to  z  =  1/2  in  the  y  =  0  plane.  The  incident 
wave  is  <£- polarized,  so  that  the  short  circuit  magnetic  field  on  the  surface  y  =  0+ 
is 

H.c(y  =  0)  =  2§-  sin  0a,e*z™9  (76) 

independent  of  (j>  for  0  <  <j>  <  tt.  This  field  induces  a  magnetic  current  Kz  in 
the  slot  having  a  functional  form  identical  to  that  given  in  Eq.  (48),  with  the 
modification  given  in  Eq.  (49);  but  the  parameter  fi  is  replaced  by  fim,  where  [8] 

ftm  =  21n(4//u;)  (77) 
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The  total  scattering  cross-section  of  the  slot  is  found  to  be  identical  to  that  of  a 
thin  wire  of  radius  w/4.  The  scattered  power  is  divided  equally  between  the  two 
sides  of  the  conducting  plane;  thus  the  transmission  cross-section  is  just  half  the 
total.  Hence 


Figure  41.  Geometry  of  the  slot-scattering  problem.  The  width  of  the  slot  is 
w  and  the  incident  magnetic  field  is  in  the  ^-direction. 

It  is  implicit  in  the  above  that  for  the  wire  or  the  slot,  the  incident  wave  is  polarized 
so  as  to  yield  maximum  scattered  power.  For  the  wire  scatterer  along  the  *-axis, 
the  electric  vector  of  the  incident  wave  is  in  the  0-direction.  From  the  results  in 
the  previous  section  one  is  led  to  conclude  that  the  equivalent  transmission  area 
of  the  slot  is  approximately  1.05  l2  at  the  fundamental  (half-wave)  resonance. 
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2.3  TIME-DOMAIN  SCATTERED  FIELD. 


We  consider  the  dipole  once  again  and  find  the  scattered  field  in  the  time  domain 
when  the  illuminating  wave  is  a  pulsed  cw  signal.  In  what  follows,  we  determine 
the  field  scattered  in  the  direction  & ,  denoting  as  before  the  direction  of  incidence 
of  the  illuminating  field  by  0.  We  have 

jkr  .1/2 

Eg  =  jkZ0sinO'-1 -  /  I(z')ejkz'C08e'dz'  (79) 

4irr  J-i/2 

Substituting  for  I(z')  from  Eqs.  (52)  -  (54)  and  carrying  out  the  integration,  we 
find 


Eil  sin  O'  t~ikr  f  2  sin  kl  sin  [y  (cos  O'  +  cos  0)] 
QD(kl)  sin  0  r  1  cos  O' +  cos  0 


(80) 


2 sin  y(l  +  cos0)j  sin  y(l  +  cos#')] 

1  +  cos  O' 

2sin  [y(l  —  cos 0)]  sin  |y(l  —  cos ^)j  ) 

1  —  cos  O'  J 

which  may  be  written  in  the  equivalent  form 

Eil  sin  O'  e~jkr  j  ( 1  +  cos  0)  cos  [kl  -  f  (cos  0  +  cos  ^)] 

E 6  QD(kl)  sin  0  r  |  (1  —  cos  O')(cos  0  +  cos  O') 

(1  —  COS  0)  COS  [&/  -I-  y  (cos  0  +  COS  S')] 

(1  4-  cos  0')(cos  0  4-  cos  O') 

2  cos  |y (cos  0  —  cost?')]  ) 
sin2  O'  J 

In  order  to  account  properly  for  causality,  we  introduce  two  modifications  to  the 
frequency-domain  result  in  Eq.  (81).  First,  we  multiply  by  the  factor  exp(—jkl  cos  0/2) 
so  that  the  wave  incident  upon  the  scatterer  will  have  zero  phase  at  the  end  of 
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the  scatterer  which  is  struck  first.  Second,  we  replace  D(kl)  by  D(kl)  exp(jkl), 
where  l/D(kl)  is  the  “pole-shifted”  Mittag-Leffler  expansion  of  the  function 


kl  sin  kl 


We  have 


D(kl)  (kl)2  jkl  ^  ( jkl  +  an)2  +  (n*-)2 


wherein  a  contribution  of  order  fl-2  has  been  neglected.  This  latter  modification 
ensures  that  the  current  at  any  point  on  the  scatterer  differs  from  zero  only  after 
that  point  has  been  illuminated. 

The  quantity  in  curly  brackets  in  Eq.  (81)  is  simply  a  weighted  sum  of  exponential 
functions.  Thus  with  the  modifications  discussed  above,  we  see  that  the  time- 
domain  scattered  field  is 


E»(r,  O', «;  0)  =  ♦  /o(<)  *  <0  «(t  -  f 


+  — -cos  0,~rfc)  (83) 

2c 


+  S(t  —  -  cos#  —  ~  cosff  —  r/c)  +  w2(0,9l)  6(t -  —  —  - cos# cos#7  —  r/c) 
c  2c  c  c  2c 

-f  6(t  + cos#7  —  r/c)  —  esc2  &  6{t  — -  — -  cos  9  +  —  cos  #'  —  r/c) 

2c  c  c  2c 


+  S(t  —  -  —  ~  cos  #7  —  r/c)  1 
c  2c  j 


where 


wx{6,&)  = 

w2(9,9')  = 


1  +  cos  9 

2(1  —  cos  #')(cos  9  -f  cos#*) 
1  —  cos  9 

2(1  +  cos#')(cos#  +  cos#') 


sin  1  U(t)  (86) 
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In  order  to  illustrate  the  time-domain  scattered  field  in  a  relatively  simple  way, 
we  consider  the  “forward-scatter”  case  wherein  O'  =  7r  —  0.  The  forward- scattered 
time-domain  field  Eg  is  given  by 

loose)  (87) 

-  £fo(t~£  cos  +  csc2  ^  W)  +  F0  (t  -  ^  j 


where 


—  F0  -  cos  0^  —  Fo  —  - - cos  0^  | 

Fo(t)  =  /  fo(t)dt 

Jo 


and  the  propagation  delay  time  r/c  has  been  omitted.  The  function  in  curly 
brackets  in  Eq.  (87)  “turns  on”  at  time  t  =  0;  its  final  contribution  begins  at  time 
t  =  2 l/c  and  decays  at  a  rate  principally  determined  by  the  damping  constant  aj, 
the  characteristic  decay  time  being  l/(a\c).  Typical  values  of  ax  are  (a  few)-1, 
so  the  total  duration  of  the  function  is  a  few  to  several  transit  times  //c.  This 
behavior  is  illustrated  in  Figures  42  to  47,  wherein  we  display  the  normalized  step 
response  (the  function  in  curly  brackets  in  Eq.  (87))  as  a  function  of  normalized 
time  ct/l  for  various  values  of  0  and  fi. 

We  may  make  use  of  this  result  to  determine  under  what  conditions  a  quasi-steady 
state  analysis  of  the  problem  can  be  applied.  If  Ei(t)  consists  of  n  cycles  of  a  cw 
signal  at  frequency  /0,  then  the  incident  pulse  length  is  n//0  sec.  A  steady-state 
condition  is  reached  in  approximately  4 Ifc  sec;  thus  n  »  4//0/c,1  the  scattered 
field  is  essentially  a  cw  pulse  of  amplitude 


Ea  = 


Mm 

nD{k0l ) 


{Ar0/sin  k0l 


—  4  csc2  6  sin  ^(l  +  cos0)  sin  —  (1  —  cos0)  l 

.  2  J  2  JJ 

'If  the  dipole  (or  slot)  is  excited  near  its  half-wave  resonance,  the  condition  becomes  n  >>  2, 
a  condition  likely  to  be  met  in  practice. 


where  k0  =  2wf0/c  and  Ei(fo)  is  the  amplitude  of  the  incident  pulse.  Further¬ 
more,  the  total  scattered  energy  can  be  accurately  approximated  by  the  product 
of  the  power  scattered  at  frequency  /0  and  the  incident  pulse  width.  Thus  un¬ 
der  conditions  where  the  slot-transmitted  or  wire-scattered  power  is  largest  (i.e., 
near  the  fundamental  resonance)  a  quasi-cw  analysis  would  appear  to  be  entirely 
appropriate. 


T 


Figure  42.  Normalized  step  response  H(t)  vs.  normalized  time  r;0  =  30°, 
0  =5. 
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H(t)  vs.  normalized  time  r; 


t 
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_ _ ap  response  H(t)  vs.  normalized  time  r; 
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Figure  47.  Normalized  step  response  H(t)  vs.  normalized  time  t\0  —  90°, 

n  =  10. 

2.4  CONCLUDING  REMARKS. 

We  have  considered  the  power  scattering  cross-section  of  an  unloaded  wire  and 
the  transmission  cross-section  of  a  narrow  slot  in  a  conducting  plane  and  have 
presented  numerical  results  to  illustrate  the  behavior  of  these  quantities  as  the 
various  problem  parameters  are  varied.  We  have  shown  that  the  equivalent  trans¬ 
mission  area  of  a  slot  of  length  l  can  be  as  large  as  (approximately)  1.05  /2,  for 
frequencies  near  its  fundamental  resonance.  We  have  also  presented  analytical 
and  numerical  results  illustrating  the  impulse  response  of  the  scattered  field  and 
have  derived  the  condition  under  which  a  steady-state  analysis  of  the  problem 
would  be  appropriate. 

We  conclude  with  a  numerical  example.  Let  the  slot  under  consideration  be  1.5  cm 
long.  Then  for  illumination  at  normal  incidence  (and  with  proper  polarization: 
E-field  across  the  slot)  and  at  a  frequency  of  10  GHz,  the  equivalent  slot  area 
will  be  approximately  2.4  cm2.  If  the  incident  power  density  is  10  watts  cm'  2, 
then  24  watts  will  penetrate  the  slot;  for  a  pulse  width  of  10  ns  (corresponding  to 
100  cycles  in  the  incident  pulse),  the  total  penetrant  energy  will  be  240  nJ.  This 
energy  will  be  radiated  into  the  system  interior. 
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SECTION  3 

PULSE  PROPAGATION  IN  WAVEGUIDES 


Summary:  We  examine  the  propagation  of  transient  signals  in  waveguides.  By 
employing  a  saddle-point  method  of  inverting  the  Laplace  transform  of  the  propa¬ 
gating  signal,  we  display  the  various  phases  of  this  signal:  the  Sommerfeld  precur¬ 
sor,  the  anterior  transient  leading  to  the  main  signal  buildup,  and  the  posterior 
transient.  We  consider  single-mode  propagation  in  some  detail;  multimode  prop¬ 
agation  can  then  be  treated  by  superposition.  As  examples,  we  investigate  the 
propagation  of  short  cw  pulses  whose  center  frequency  is  above  the  dominant¬ 
mode  cutoff  frequency  and  whose  envelopes  are  initially  Gaussian  or  rectangular. 

The  propagation  of  pulsed  signals  in  waveguides  is  a  subject  of  critical  importance 
for  studying  the  potential  effects  of  microwave  weapons  on  systems.  The  subject 
is  not  a  new  one,  especially  in  its  general  form  wherein  the  mathematical  problem 
is  that  of  treating  propagation  in  a  dispersive  medium.  A  waveguide  constitutes 
a  special  case  of  the  more  general  problem. 

We  shall  address  the  problem  with  a  view  to  demonstrating  the  principal  features 
of  the  dispersed  signal.  We  draw  most  heavily  on  the  work  of  Sommerfeld  [9] 
and  Haskell  and  Case  [10];  Sommerfeld  investigated  the  percursor  signal  which 
propagates  in  a  dispersive  medium  and  which  is  the  first  signal  component  ob¬ 
served  at  early  retarded  times,  while  Haskell  and  Case  considered  transient  signal 
propagation  in  a  cold  plasma.  The  latter  problem  has  a  mathematical  form  which 
corresponds  to  the  waveguide  propagation  problem. 

The  propagation  of  a  cw  signal  of  frequency  w  in  a  given  waveguide  mode  gives 
rise  to  the  simple  relation 

f(z,t )  =  Acos(ud  —  -\Ju)2  —  oj]  +  <j>)  (90) 

where  f(z ,  t )  is  the  signal  at  a  point  z  in  the  guide  and  A  and  <j>  are  the  amplitude 
and  phase  of  the  signal  at  z  =  0.  It  is  assumed  in  (90)  that  u  >  u)c,  the  cutoff 
frequency  of  the  mode  under  consideration.  In  the  event  that  w  <  wc,  we  have 

f(z,t)  =  Acos(ut  +  ^)exp  —  u>2^  (91) 
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and  the  signal  does  not  propagate.  The  phase  and  group  velocities  of  a  monochro¬ 
matic  signal  are  defined 

Up  =  - - -  (92) 

yj\  -  (Wc/w)2 

vg  =  C\J\  -  ( u>c/u ;)2  (93) 

where  c  denotes  the  speed  of  light.  The  fact  that  these  velocities  are  not  indepen¬ 
dent  of  the  signal  frequency  indicates  that  the  guide  is  dispersive:  a  pulsed  signal 
will  be  spread,  or  dispersed,  since  its  higher-frequency  components  travel  with  a 
greater  group  velocity  than  the  lower-frequency  components. 

The  signal  dispersion  can  be  illustrated  by  considering  the  propagation  of  a  Gaus¬ 
sian  pulse  of  center  frequency  u;0  (u;0  >  u>c)  and  width  t0: 

/(0,  t)  =  Ae~t2/2t°  cosu>0t  (94) 

The  signal  at  a  point  z  >  0  in  the  guide  can  be  written  as  a  Fourier  integral: 


/(M)  =  ^  h(. ;'w)eJ(u't 


f  x/w2—  u/?  . 


where 


Fo(ju)  =  AtoyJ |  [e-("-^)2‘o/2  +  e-(w+wo)%’/2j 


By  expanding  the  function 


g(ju)  =j(ut-^y 


'u2  —  u>2 


near  u>  =  lo0  as 


g(joj)  =  ju0(t  -  z/vv o)  +j(u>  -  u0)(t  -  z/vgo )  +  -j(u>  -  u>0)2p  (98) 

where  Vpo  and  Vgo  are  as  given  in  (92)  and  (93)  with  u  =  u>o, 


C  (u0!-U>?)W 
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and  similarly  for  u  =  —u>0  and  evaluating  the  integral  in  (95),  we  find 


/(M) 


=  A 


(t  -  Z/Vg0)2 

2*8(1  +  p7<$) 


(100) 


cos 


u>o{t  -  z/Vpo)  +  -  tan 


p(t  -  z/vg q)2 
2*8(1  +  P2/*8). 


We  remark  that  if  the  term  in  (98)  involving  the  parameter  p  were  neglected,  we 
would  obtain 

/(*,  o  s  Ae-it-z/v’°W2t l  cos  _  z/Vpo)]  (101) 

indicating  that  the  pulse  envelope  travels  at  the  group  velocity  vg0  and  the  points 
of  constant  phase  travel  at  the  phase  velocity  VpQ.  Keeping  the  term  involving  p, 
we  see  that  while  the  center  of  the  pulse  envelope  travels  at  speed  Vgo,  the  width  of 
the  envelope  increases  as  it  propagates,  and  its  amplitude  decreases.  Furthermore, 
the  phase  changes  in  a  complicated  way  as  the  pulse  propagates. 

We  note  that  the  governing  parameter  in  this  example  is 


P 

tl 


(102) 


from  which  it  is  evident  that  pulse  distortion  will  occur  if  any  of  the  following 
conditions  are  met: 


•  the  propagation  distance  (measured  in  units  of  ct0)  becomes  large; 

•  the  pulse  width  (measured  in  units  of  u;"1)  becomes  small;  or 

•  the  center  frequency  of  the  pulse  becomes  close  the  the  cutoff  frequency. 

We  show  in  Figure  48  the  variation  in  the  pulse  envelope  as  the  pulse  propagates 
as  a  function  of  normalized  retarded  time  ( t  —  z/vgo)/to  for  different  values  of 
normalized  propagation  distance  p/t\. 

Having  demonstrated  the  phenomenon  with  this  simple  example,  we  turn  in  sub¬ 
sequent  sections  of  this  chapter  to  a  more  general  treatment  of  pulse  dispersion 
in  waveguides.  We  examine  single-mode  transient  propagation  and  then  discuss 
a  simple  generalization  to  the  multimode  case. 
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Figure  48.  Gaussian  pulse  envelope  vs.  normalized  time  ( t  —  z/vgo)/t0  for 
various  values  of  normalized  distance  p/t% . 

3.1  FORMULATION:  SINGLE-MODE  PROPAGATION. 

The  calculation  of  the  dispersed  signal  associated  with  the  propagation  of  a  single 
waveguide  mode  of  radian  cutoff  '"equency  u>c  reduces  to  the  problem  of  evaluating 
the  integral 

/(M)  =  ST  (io3) 

where  T  is  the  Bromwich  contour  (a  —j oo  to  cr+joo)  in  the  complex  s-plane  and 
where  Fo(s)  denotes  the  Laplace  transform  of  /(0,  t): 

Fo(s)=  /°°  /(0,  (104) 

Jo 

We  shall  assume  in  the  following  that  the  only  singularities  of  Fo{s)  are  poles  in 
the  left  half-plane.  The  points  s  =  ±  jwc  are  joined  by  a  branch  cut  so  that  the 
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5-plane  appears  as  shown  in  Figure  49.  The  cut  is  chosen  so  that  on  the  upper 
Riemann  sheet, 


(s2  +  w2)1/2  =  JRiRi exp \j(fa  +  fa)  12]  (— 7r/2  <  fa,  fa<  37r/2)  (105) 


We  shall  evaluate  f(z,t)  approximately,  using  an  appropriate  set  of  analytical 
techniques.  We  shall  show  that  in  many  cases  of  interest,  the  dispersed  signal 
comprises  three  phases:  first,  the  “Sommerfeld  precursor”  which  appears  at  early 
times  ( t  slightly  greater  than  z/c);  next,  an  anterior  transient  phase  which  precedes 
the  main  signal  buildup  at  times  near  z/u3,  where  vg  is  the  group  velocity  near 
the  dominant  frequency  of  /(0,t);  and  finally,  a  posterior  transient  phase  at  late 


Figure  49.  The  complex  s-plane.  The  x’s  indicate  poles  of  the  function  F0(s). 


3.1.1  The  Sommerfeld  Precursor. 

We  begin  by  deforming  the  Bromwich  contour  into  a  semicircle  of  (large)  radius 
R  in  the  right  half-plane  and  adding  its  image  in  the  left  half-plane  as  shown  in 
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Figure  50,  so  that 

f(z,  t)  =  ~  [  Fo(s)eat~i\^  ds  (106) 


The  contribution  from  the  closure  in  the  left  half-plane  is  negligible  for  t  >  z/c; 
Cr  denotes  the  circular  integration  path.  On  the  path,  we  assume  that 


(->1) 

sn 

(107) 

where  A  is  a  constant,  and  write 

2 

st  -  -Js*  +  u;2  =  a(<  -  z/c)  -  +  •  •  • 

(108) 

Thus 

/(*,()  =  A-.l 

v  ’  ’  27 r;  Vch 

(109) 

where  r  =  t  — 

z/c  is  the  retarded  time  and 

III 

(110) 

Now  the  integral  in  (109)  can  be  evaluated  by  defining  a  new  variable  <f>  through 


(111) 


and  identifying  yjifr  as  the  radius  of  the  circular  integration  path.  Then 


/(M) 


.  /  \  (n-l)/2 

—jn~l  (  -  )  P  d<j> 

2 7T  \^  /  J  —  7T 


(112) 


=  A 


2 (ct  -  z) 


.  (n-l)/2 


dfl—l 
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Figure  50.  Deformation  of  T  into  a  contour  suitable  for  evaluating  the 
Sommerfeld  precursor.  CR  denotes  the  circular  contour. 

where  Jn-i(-)  denotes  a  Bessel  function.  Evidently,  if  Fo{s)  can  be  expanded  in  a 
series  of  the  form  ^ 

^o(s)  =  53"r  (S-MX>)  (113) 


then  the  precursor  signal  can  be  expressed  as 


00  \2(ct  —  r- 

/(«.*)  ^ 

1=1 


U>c\  ~  [t  -  - 


(114) 


at  early  times. 


To  illustrate  the  application  of  this  result,  we  consider  a  specific  example.  Let 


f(0,t)  =  a0  sin u>0t  U(t) 


so  that 


M>)  = 


s3  +  u>0 


(115) 


(116) 


(VqCIq  UJqCLq 


Thus  the  Sommerfeld  precursor  is 


''•'■'-Ha's-oP'H? 


c  \  z 


(117) 


which  involves  the  frequency  t oQ  only  in  the  amplitude  factor.  The  oscillatory 
character  of  this  function  depends  only  upon  the  cutoff  frequency  u)c.  We  shall 
return  to  consider  the  Sommerfeld  precursor  after  we  examine  the  next  phase  of 
the  dispersed  signal,  the  anterior  transient. 


3.1.2  The  Anterior  and  Posterior  Transients. 

To  investigate  the  behavior  of  the  dispersed  signal  after  the  earliest  times,  we 
employ  a  saddle  point  integration  technique.  Define 


g(s)  =  st - v/s2  +  u/2 

c  v 


This  function  has  two  saddle  points,  located  at 


S±  =  ±ju,p  =  ±ju]c 


ict/z) 

f(ct/z)2  —  1 


(118) 


(119) 


We  expand  y(s)  about  the  saddle  points  in  Taylor  series,  obtaining 


g{s)  =  g(±jujsp)  +  ^  (-s  T  ;'w,P)2  + 

*± 


(120) 
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where 


#9 


=F  jz 


ds 2 


a=±ju,p 


U)cC 


(121) 


Now  deform  the  original  contour  T  into  the  steepest-descent  path  SDP,  as  shown 
in  Figure  51.  We  have  shown  in  the  Figure  the  fact  that  the  deformation  of  T 
into  SDP  can  pick  up  some  of  the  poles  F0(s)  lying  in  the  left  half  of  the  s-plane. 


Thus 


/(z,<)  =  - ^  /  F0(s)e9{a)ds  +  £  Res  £>(*)< 

2lCJ  JSDP  n 


,9(a) 


(122) 


where  sn  denotes  a  pole  crossed  by  the  deformation  of  the  contour.  By  virtue  of 
the  fact  that  the  poles  crossed  will  occur  in  the  complex-conjugate  pairs  (sn  = 
-an  ±  jun)  and  that  f(0,t)  is  a  real  function,  we  have  for  a  typical  pole-pair 


contribution 


fpn(t ,  z)  =  2^„e-ff"t+Q"2  cos(u :nt  -  Pnz  +  <Pn)U(t  -  <„)  (123) 


where 

Res  F0(s) |_„B+J Wn  =  Rpne (124) 

-  Qn  +  jPn  =  \/^c  ~  +  j<*n)2  (125) 

t  =  *  Un  ~ g"  (126) 

"  C  \J (Wn  —  On)2  — 

When  the  saddle  points  are  not  close  to  poles  of  Fq[s),  the  contributions  to  /(z,t) 
from  the  neighborhoods  of  the  two  saddle  points  are  easily  found.  We  obtain  for 
/ap(z,<),  the  contribution  to  f{z,t)  from  integration  along  SDP, 


/*p(M)  = 


(127) 


cos 


-  r+  -  +  it 
4 
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Figure  51.  The  steepest-descent  contour  in  the  5-plane.  Poles  in  the  region 
bounded  by  the  dashed  lines  cannot  be  crossed  by  SDP. 


where  we  have  written 


Fo{ju,p) 


Fo  (jw,p) 


(128) 


The  above  result  applies  when  the  saddle  points  are  neither  too  close  to  poles 
of  F0(s)  which  might  lie  close  to  the  ju;-axis  nor  too  close  to  the  branch  points 
at  s  =  ±ju>c,  which  are  the  points  approached  by  the  saddle  points  at  very  late 
times.  Before  proceeding  to  the  investigation  of  these  situations,  however,  we 
observe  that  if  F0(s)  as  given  in  our  example  for  the  Sommerfeld  precursor  is 
applied  to  the  calculation  of  f,p(z,t)  at  early  times,  we  find 


r  (  ( z  V  (  2c  V/2  (ct 

Jsp{Z,t)  -*•  AO  I  —  ]  -  )  l  -  I  ~ 

\U>C )  \ctj  \TTUJ-zJ  \Z 


1/4 


(129) 


cos 
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when  ct/z  is  only  slightly  greater  than  unity.  One  immediately  observes  that 
the  result  in  (129)  is  identical  to  that  obtained  if  one  uses  (117)  with  the  usual 
asymptotic  approximation  for  Thus  the  early-time  saddle  point  integral 

expression  for  f(z,t)  matches  the  Sommerfeld  precursor  expression  when  ct/z  — 
1  <<  1  and 

'/5v\/f^>>1  (130) 

The  required  conditions  can  be  simultaneously  met  if  z  is  at  least  a  few  wave¬ 
lengths  at  the  cutoff  frequency.  We  shall  assume  that  2  is  always  at  least  this 
large,  and  we  identify  the  transition  between  the  Sommerfeld  precursor  and  the 
early-time  saddle-point  expression  as  the  onset  of  the  anterior  transient. 

When  the  saddle  points  pass  close  to  a  complex-conjugate  pole  pair  the  expression 
for  fsp(z,t)  given  in  (127)  must  be  modified.  We  remark  parenthetically  that  if 
the  poles  are  located  at  s  =  —  av  ±  jup  where  ap  is  small  and  u>p  >  u>c,  the  saddle 
points  and  the  poles  are  closest  when  u>sp  =  u>p,  or  when 

2  =  ct  yj  1  (u;c/u>p)2  (131) 

= 


where  vg(ujp)  is  simply  the  group  velocity  for  a  signal  at  frequency  u>p. 
When  a  pole  is  close  to  a  saddle  point,  we  write  for  Fq{s) 


F0(s)  = 


Res  F0(sp) 


S  «S*i 


+  Go(s) 


(132) 


where  Go(s)  is  analytic  at  s  =  sp.  Then  the  contribution  to  f(z,t)  from  the 
integration  over  SDP  can  be  written 


1 

2ttj 


L 


SDP 


F0{s)eais)ds 


2  nj 


Res  Ft 


o(sP)Js 


Pp(») . 


ds 


SDP 


(133) 
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+  7^-1  G0(s)e^ds 

2irj  Jsdp 

The  second  integral  on  the  right-hand  side  of  (133)  yields  a  contribution  of  the 
form  given  in  (127)  and  (128)  with  Fo  replaced  by  G0-  The  first  integral  can  be 
expressed  in  terms  of  the  complementary  error  function;  and  we  find  that  the 
contribution  to  the  integral  over  SDP  made  by  poles  near  the  saddle  points  before 
SDP  crosses  over  the  poles  is 

/,P(M)  =  (134) 


jRes  Fo(<)ei(“‘~‘v^)e (o 


^Itss  Fc,(sr)e  J'/4^(“w -J4,)] 


where  s±  =  —<rp  ±  ju>p  and 


2  13/2 


(135) 


After  SDP  has  passed  over  the  poles,  we  include  the  pole  contributions  to  obtain 
/(*,<)  =  ^Res  ^(s+)e;(^_fv^)  •  (136) 

[2  -  erfc  [eW4y!  (-;<  -  Wjp)]  j 

+  iResF(s;)e-J(-‘^v^F). 

{2  -  erfc  (js;  -  u>ip)]  } 

We  observe  that  f(z,t)  as  given  in  (134)  and  (136)  above  is  continuous  as  SDP 
picks  up  the  pole  contributions;  this  is  a  consequence  of  the  symmetry  property 
of  the  complementary  error  function  (11]: 


erfc(— iz)  —  2  —  erfc(iz) 


(137) 
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At  very  late  times,  the  saddle  points  approach  the  branch  points.  It  is,  however, 
not  necessary  to  modify  the  result  in  (127)  in  the  limiting  case,  since  the  magnitude 
of  g"(s)  becomes  very  large  as  the  branch  points  are  approached.  Thus 

'2u ;cc\1/2  /z\3/2 


f,P(z,  t)  ~  (^)  |^o(M)|  cos  (u >ct  +  j  (138) 


as  t  — >  oo. 


We  now  consider  our  example  once  again,  where  /( 0,  t)  is  a  sine  wave  of  frequency 
u;o  turned  on  at  t  =  0.  The  precursor  signal  which  appears  at  the  earliest  times  is 

=  (139) 

The  anterior  transient  signal  before  the  saddle  points  approach  the  poles  at  ±jw0 
is 


..  /2wc\1/2  [7cA2  V'U/cAV  2  2\  2 

f(z it)  —  aouJo  —  ^  \T/  \  c  ~  wo)  +  uo 


-l 


(140) 


cos  I  wc 


im7'-*) 


When  the  saddle  points  approach  and  pass  through  the  poles,  we  have  for  wsp  >  w4 
(i.e.,  before  the  crossing  occurs) 


/(z,()  =  0  V^F^)erfc  [e>-/-  y|(az„  -  Wo) 

-  c-J(u'ol_‘ v/“°~“'0erfc  [e-»-/<y|Kp  -  „„)]  | 


(141) 


After  the  poles  are  crossed, 


/(*,<)  =  [2-erfc(e'''’y|K-u,w))]  (142) 

[2  -  erfc  (e-W^(uo  -  „„))]} 


—  e 
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Finally,  after  the  saddle  points  have  moved  away  from  the  poles,  we  have 
cos  ~  1  +  j)  +  oo sin  (u0t  -  ^\/wo  j 


u>0t  -  ~\jul~u2c 
\  c 


which  comprises  the  posterior  transient  and  the  steady-state  signal  at  frequency 
u>o- 

Near  the  time  of  the  main  signal  buildup,  the  function  in  (141)  can  be 

written  in  the  form 


f{z,t)  =  -sinw0(i  -  z/vpo) 

1  ~  (_)yn-H  r 


(144) 


1  ~  -)n^"+1  .  ,  .  .  Tilt  7T 

+  ~7=  Y,  to  I  Sm  ~  2  VT*>)  +  +  J 

V7rn^(2n  +  1)n!  2  4 


wherein  Upo  denotes  the  phase  velocity  at  frequency  w0: 

Vpo  =  c(l  -u2Ju2)~1'2 


(ljcz\1/2  (ct\2 

q=(^J  [Vlj  ~ 


ct  u>o  /  ( cV 


1  - +  - 


(145) 


(146) 


One  notes  that  q  passes  through  zero  when  t  =  z/vg 0,  where  Vgo  is  the  group 
velocity  at  frequency  u>o 

vgo  =  c(l  -u2clul)  ,2  (147) 

If  we  write 

/(z,  t)  =  A(z,  t)  sin  [u^0(<  -  zfapo)  +  <f>(z,  0]  (148) 

where  A{z,t)  is  an  amplitude  and  <j>(z,t)  a  phase  function,  we  can  readily  show 
that  the  amplitude  function  A(z,  t )  is  given  by 


^(*,o  =  [G  +  5,)2+'!?2’ 


(149) 


62 


where 

C  _  J_  Y'  (~)Vn+1  COS  /nr 
2  v/ir  (2n  +  l)n!  sin  \  2  4/ 

When  g  =  0,  we  have 


(150) 


A(z,t) 

whence  the  rate  of 
buildup  is  clearly 


<i5i) 

rise  of  the  signal  amplitude  near  the  time  of  the  main  signal 


dA 

dt 


Uc 

2 


(152) 


t  =  Z/Vg 0 

The  decrease  in  the  rate  of  rise  as  kcz  increases  or  as  u0  approaches  u;c  is  apparent. 

We  show  in  Figure  52  a  plot  of  /(z,t)  as  a  function  of  normalized  retarded  time 
rr  =  u>0 (t  —  z/c)  for  the  case  u)o/u>c  =  1.5,  kcz  =  60tt.  The  saddle  point  crosses  the 
pole  at  rr  =  96.6  for  this  parameter  choice.  One  observes  that  the  precursor  and 
anterior  transient  signals  build  continuously  in  amplitude  from  rr  =  0  and  that 
the  sum  of  the  pole  contribution  and  the  posterior  transient  produce  a  slow  “beat” 
in  the  envelope  of  the  signal  after  the  pole  has  been  crossed.  The  amplitude  of  the 
signal  gradually  approaches  unity  at  late  times,  the  beat  amplitude  decreasing  as 

r3/2. 


3.2  A  PULSE-PROPAGATION  EXAMPLE. 

We  now  consider  the  application  of  the  foregoing  analysis  to  a  practically  inter¬ 
esting  problem  of  pulse  propagation.  Specifically,  we  assume  that 

/(0,  t)  =  a0smu>0t  [t/(0  —  U{t  —  t0)]  (153) 

where  Woto/71-  is  equal  to  an  even  integer. 
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Turn-On  Sine  Wave 


Figure  52.  Turn-on  sine  wave  of  frequency  u>0  as  a  function  of  normalized 
retarded  time  u}0(t  —  z/c)  for  u>0/uc=l. 5;  normalized 
propagation  distance  kcz  =  607r. 


By  virtue  of  the  fact  that  this  signal  is  simply  the  difference  of  two  turn-on  sine 
waves  separated  by  an  interval  to,  we  have  immediately 


/(•M)  =  fo(z,t)  -  f0(z,t  -  t0)  (154) 

where  /0(z,  t)  now  denotes  the  function  given  in  (139)  to  (143)  above. 

We  show  in  Figures  53  -  58  the  (pulse)  function  f(z,t )  as  a  function  of  normalized 
retarded  time  u>0(t  —  z/c)  fo;  u;0/wc  =1.5,  u>0t  =  8ir  and  various  values  of  kcz.  The 
results  show  the  dispersion  of  the  pulse  as  a  function  of  propagation  distance. 
We  note  that  when  kcz  is  not  large,  the  peak  value  of  the  propagated  signal  can 
exceed  its  value  at  z  =  0.  This  is  related  to  the  “beat”  phenomenon  mentioned 
earlier  in  connection  with  the  turn-on  sine  wave. 
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Sine-Wave  Pulse 
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Figure  53.  Sine- wave  pulse  vs.  normalized  retarded  time  u>o(t  —  z/c)  for  kcz  —  0 
C Jq^O  ~  8tT- 


Sine-Wave  Pulse 


Figure  54.  Sine-wave  pulse  vs.  normalized  retarded  time  u>o {t  —  z/c)  for 
kcz  =  57r;  u;0/o  =  S7T . 
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Figure  57.  Sine-wave  pulse  vs.  normalized  retarded  time  u>0{t  —  zjc)  for 
kcz  =  40?r;  LOOt0  —  8 ir. 

Sine-Wave  Pulse 


Figure  58.  Sine- wave  pulse  vs.  normalized  retarded  time  u0(t  —  z/c)  for 
kcz  =  607r;  u>0 t0  =  87r. 
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3.3  MULTIMODE  PROPAGATION. 


The  transverse  electric  field  inside  a  conducting  waveguide  of  general  cross-section 
can  be  expressed  in  terms  of  the  individual  modal  field  distributions  emn(£,7)  in 
the  form 

Et({,  t)  =  amnemn(£,  r/)fmn(z,  t )  (155) 

m,n 

where  amn  is  the  amplitude  of  the  (m,n)th  mode  and 

/«„(*,  t)  =  ~  jf  Fo(s)  exp  (st  -  ~cyJ S2  +  u2cmr}j  ds  (156) 

The  frequency  u>cmn  is  the  cutoff  frequency  of  the  (m,  n)th  mode  and  £  and  rj  are 
the  transverse  coordinates.  It  is  assumed  that  the  time  dependence  /( 0,  t)  is  the 
same  for  each  mode.  Each  of  the  integrals  (156)  can  be  treated  as  in  the  foregoing 
and  the  results  superposed  to  yield  the  total  transverse  electric  field  at  any  point 
in  the  guide. 

We  consider  a  simple  example  to  illustrate  the  multimode  propagation  problem. 
Let  the  waveguide  be  rectangular  in  cross-section  such  that  the  interior  is  defined 
by  0  <  x  <  a,  0  <  y  <  b,  and  let  the  electric  field  at  z  =  0  be  E0f{t)ay.  Then  in 
the  guide,  the  electric  field  Ey  for  any  position  z  >  0  is 


Fy{x,y,z;t) 


Z 


\Eq  .  mrx 
- sin - 


n  w  a 


(odd) 


(157) 


^7  Jr  F0{s )  exp  (st  -  ds 

where  uicn  =  nwc/a.  Consider  the  field  across  the  center  of  the  guide  (x  =  a/2) 
when  f(t)  is  a  sinusoidal  pulse  with  a  Gaussian  envelope.  Then  from  results 
obtained  earlier  in  this  section,  we  obtain 


fa  \  "m“ 

\2  ’  Z'  /  =  ^ 


n«l 

(odd) 


4  E0  ,  nir 

- sin  — 

nir  2 


(158) 


exp 


~(<  ~  Z/Vgrl )2 

2tl(l+P2Jt4o). 


cos  [w0  ( t  -  z/vp n)  +  <t>n{z,  t)] 
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where  1 0  is  a  measure  of  the  pulse  width  at  t  =  0  and 


Pn 


v 


gn 


V 


pn 


<t>n{z,t) 


c^M 

<?/vgn 

—  tan-1  —  —  P"  ~  z/vgn) 

2  tl  2t$(l  +p»J?0) 


(159) 


and  Umax  is  the  largest  value  of  n  such  that  o>cn  <  ojq. 

We  have  computed  EyJE0  as  a  function  of  normalized  retarded  time  u;0(t  —  z/c) 
for  various  values  of  normalized  2-position,  for  the  case  u0t0  =  4tt,  u?0  =  5.5tt c/a. 
The  normalized  distance  is  taken  to  be  kcXz  =  u>clz/c.  Three  propagating  modes 
(n  =  1,  3,  5)  are  excited  in  the  waveguide.  Numerical  results  are  presented  in 
Figures  59  to  64.  Multimode  Gauss  ion  Pulse 


Figure  59.  Mode  sum  Ey/ E0  vs.  normalized  retarded  time  w0(t  —  z/c)  for 
Gaussian-envelope  pulse.  Normalized  propagation  distance 
kcXz  =  0;  ujotQ  =  4?r. 
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Figure  60.  Mode  sum  Ey/E0  vs.  normalized  retarded  time  u >0{t  —  zjc)  for 
Gaussian-envelope  pulse.  Normalized  propagation  distance 
kc\Z  —  u?o^o  = 

Multimode  Gaussian  Pulse 


Figure  61.  Mode  sum  Ey/E0  vs.  normalized  retarded  time  u>0{t  —  z/c)  for 
Gaussian-envelope  pulse.  Normalized  propagation  distance 
kcXz  =  5 ff;  u>0t0  =  4?r. 


Multimode  Gaussion  Pulse 


Figure  62.  Mode  sum  Ey/Eo  vs.  normalized  retarded  time  u>o [t  ~  z/c)  for 
Gaussianenvelope  pulse.  Normalized  propagation  distance 
kciz  =  107r;  w0t0  =  47T. 

Multimode  Gaussian  Pulse 


Figure  63.  Mode  sum  Ev/Eo  vs.  normalized  retarded  time  u>o {t  —  z/c )  for 
Gaussian-envelope  pulse.  Normalized  propagation  distance 
kciz  =  20tt;  utoto  =  4?r. 
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Multimode  Gaussian  Pulse 


Figure  64.  Mode  sum  Ey/E0  vs.  normalized  retarded  time  u0(t  —  z/c)  for 
Gaussian-envelope  pulse.  Normalized  propagation  distance 
kciz  =  307r;  u>0t0  =  47T. 


Figure  59  shows  the  sum  of  the  three  modes  at  z  =  0.  The  peak  amplitude  exceeds 
unity  because  the  non-propagating  modes  have  not  been  included  in  the  sum.  In 
Figures  60  to  64  we  show  the  signal  at  different  positions  in  the  guide.  One  will 
note  that  as  the  distance  kcXz  increases  the  TE50  mode  contribution  separates 
first,  followed  by  the  TE30  contribution.  When  kclz  =  5n,  the  TE50  contribution 
is  already  very  highly  dispersed,  while  the  TE30  contribution  is  just  beginning  to 
appear.  The  development  and  separation  of  this  contribution  is  evident  in  the 
remainder  of  the  plots.  It  is  interesting  to  note  that  after  the  higher-order  mode 
contributions  have  been  separated,  the  dominant-mode  contribution  exceeds  the 
three-mode  sum  at  z  =  0.  This  is  a  consequence  of  the  fact  that  the  TEjo  mode 
amplitude  coefficient  is  4/n  =  1.27;  when  this  modal  contribution  is  separated  from 
the  higher-order  modes  by  the  dispersion,  the  total  signal  amplitude  displays  the 
peak  value  of  the  TE10  mode  alone. 

Figures  65  -  70  display  the  three-mode  sum  for  the  case  u>0t0  ~  27T.  The  increased 
dispersion  associated  with  the  broader  bandwith  of  this  pulse  is  apparent. 
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Multimode  Gaussian  Pulse 


Figure  65.  Mode  sum  Ey/E0  vs.  normalized  retarded  time  u>0{t  —  z/c)  for 
Gaussian-envelope  pulse.  Normalized  propagation  distance 
kciz  =  0;  u0t0  =  2n. 

Multimode  Gaussian  Pulse 


Figure  66.  Mode  sum  Ev/Eo  vs.  normalized  retarded  time  u>0(t  —  z/c)  for 
Gaussian-envelope  pulse.  Normalized  propagation  distance 
kc\z  =  27r;  u)0tQ  =  2n. 
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Multimode  Gaussian  Pulse 


Figure  67.  Mode  sum  Ey/E0  vs.  normalized  retarded  time  u0(t  —  z/c)  for 
Gaussian-envelope  pulse.  Normalized  propagation  distance 
kclz  —  5tt;  u>0t0  =  2 7r. 

Multimode  Gaussian  Pulse 


kcXz  =  IOt;  uJoto  =  2ir. 
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Multimode  Gaussian  Pulse 


Nornaiized  Retar dad  Time 


Figure  69.  Mode  sum  Ey/E0  vs.  normalized  retarded  time  u0(t  —  z/c)  for 
Gaussian-envelope  pulse.  Normalized  propagation  distance 
fccl  z  =  20tt;  u)0tQ  =  27r. 


Multimode  Gaussian  Pulse 


Norealized  Retarded  Ttee 


Figure  70.  Mode  sum  Ev/E0  vs.  normalized  retarded  time  u0(t  -  z/c)  for 
Gaussian-envelope  pulse.  Normalized  propagation  distance 
kc\z  =  307r;  uqIq  =  27t. 
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3.4  CONCLUDING  REMARKS. 

In  this  chapter  we  have  considered  the  approximate  analytical  evaluation  of  tran¬ 
sient  signals  propagating  in  waveguides.  Both  Gaussian  and  rectangular  pulses 
have  been  used  as  examples.  The  multimode  problem  was  also  discussed  and 
illustrated  by  example. 

The  propagation  of  a  pulse  whose  center  frequency  is  not  much  greater  than  the 
fundamental  mode  cutoff  frequency  is  characterized  by  substantial  dispersion.  The 
peak  amplitude  and  the  rate  of  rise  of  the  pulse  decrease  as  it  propagates,  while  the 
apparent  pulse  length  increases.  A  pulse  whose  center  frequency  greatly  exceeds 
the  cutoff  frequency  propagates  with  little  distortion;  but  several  modes  would 
typically  exist  simultaneously  in  such  a  situation.  It  was  shown  in  an  example 
that  the  dominant  mode  propagated  almost  undistorted  and  that  the  appearance 
of  the  total  signal  was  essentially  that  of  the  dominant  mode  followed  by  a  “tail” 
comprising  the  higher-order  mode  contributions.  Thus  when  the  center  frequency 
of  the  pulse  is  well  above  the  dominant-mode  cutoff  frequency,  the  dominant¬ 
mode  contribution  to  the  total  field  will  be  the  first  to  arrive  at  a  given  point 
and  will  likely  be  the  largest-amplitude  portion  of  the  total  signal.  Multimode 
propagation  per  se  would  not  appear  to  be  an  issue  in  microwave  interactions  with 
systems;  the  important  point  would  seem  to  be  simply  that  a  signal  well  above 
the  dominant-mode  cutoff  frequency  can  propagate  essentially  undistorted  in  the 
system  waveguides. 


SECTION  4 

POWER  COLLECTED  BY  AN  ANTENNA  OVER  A  LOSSY 

ROUGH  SURFACE 


Summary:  We  develop  an  expression  for  the  power  received  by  an  antenna  over  a 
lossy  rough  surface.  We  perform  representative  calculations  for  the  case  where  the 
illumination  is  randomly  polarized  and  normally  incident,  and  where  the  rough- 
surface  height  is  a  two-dimensional  stationary  random  process  with  isotropic  cor¬ 
relation. 

A  receiving  antenna  situated  over  a  rough  surface  which  is  illuminated  by  a  re¬ 
mote  transmitter  will  receive  power  both  directly  from  the  illuminating  field  and 
indirectly  from  the  field  reflected  by  the  surface.  The  reflected  field  will  comprise 
both  coherent  and  incoherent  contributions,  the  latter  resulting  solely  from  the 
roughness  of  the  surface.  The  coherent  contribution  arises  from  reflection  from  the 
“average”  surface.  If  the  surface  height  variations  are  of  the  order  of  the  incident 
wavelength  or  larger,  the  coherent  contribution  is  very  small  and  the  reflected  field 
is  dominated  by  the  incoherently  scattered  contribution. 

We  consider  the  power  received  from  the  incoherently  scattered  field.  Although 
rigorous  theories  of  scattering  from  very  rough  surfaces  do  not  exist,  an  approxi¬ 
mate  theory  based  on  first-order  perturbation  methods  does  exist  and  can  be  used 
to  address  the  problem.  The  approximate  theory  is  strictly  valid  only  when  the 
surface  height  variations  are  small  in  comparison  to  the  wavelength,  but  it  pro¬ 
vides  a  phenomenological  framework  within  which  the  more  general  case  can  be 
studied.  We  employ  the  first-order  perturbation  theory  in  this  section  and  remark 
as  appropriate  on  the  use,  for  example,  of  experimentally  derived  scattering  data 
in  the  theoretical  framework. 

In  the  following  we  review  the  rough-surface  scattering  problem  as  it  applies  to 
a  rough  interface  between  free  space  and  a  lossy  medium.  Then  we  consider  the 
power  collected  by  an  antenna  over  such  an  interface  when  the  illumination  is 
normally  incident  from  above.  We  end  the  section  v/ith  an  example  calculation. 
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4.1  ELECTROMAGNETIC  SCATTERING  FROM  A  LOSSY 
ROUGH  SURFACE. 

Consider  a  monochromatic  plane  electromagnetic  wave  incident  at  an  angle  0-  on 
a  rough  surface,  as  shown  in  Figure  71.  The  incident  wave  carries  power  densities 

•  Sfnc  in  horizontal  polarization 

•  S?nc  in  vertical  polarization 

The  power  densities  in  each  polarization  which  are  incoherently  scattered  by  an 
area  A  A  on  the  rough  surface  in  the  direction  9S,  <j>3  are  given  by 


=  AA 

4  7rr2 


&hh  °hv 


&vh 


(160) 


where  crhh, . . .  denote  the  scattering  cross-sections  per  unit  area  of  the  rough  sur¬ 
face  and  R  is  the  distance  from  the  area  AA.  Each  of  these  cross-sections  depends 
on  the  angles  0,,  9„,  and  <j)3\  the  (complex)  relative  permittivity  eT  of  the  scattering 
medium;  the  frequency;  and  the  statistical  properties  of  the  rough-surface  height. 
Specifically,  the  first-order  perturbation  theory  yields  [12] 

<7  =  Airk4  cos2  0,  cos2  0aW(fcsin  9,  cos  (j),  —  k  sin  0,,  k  sin  0,  sin  |  a  |2  (161) 

where  k  is  the  propagation  constant  of  a  plane  wave  in  free  space  and  where 


<*hk  = 


Othv  — 


&vh 


_ ~(Cr  ~  1) COS  (ft, _ 

^cos  0,  +  \JtT  —  sin2  0,  j  ^cos  9,  -f  \Jcr  —  sin2  0,^ 

_ (er  ~  l)sin<J>3\/eT  -  si n20, _ 

^er  cos  0,  +  \]t.T  —  sin2  0,^  ^cos  0,  +  \JtT  —  sin2  0,^ 

_ ~(cr  -  1)  sin  <t>,\jtT  -  sin2  0, _ 

^cos  0,  +  \jtT  —  sin2  0,j  (ytT  cos  93  +  \JeT  —  sin2  0,y 


(162) 


(163) 


(164) 
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area  AA 


Figure  71.  Geometry  of  the  rough-surface  scattering  problem. 

(er  —  1)  ^er  sin  0*  sin  9S  —  cos  < pa\JtT  —  sin2  Qi\Jtr  —  sin2 


(165) 


^er  cos  0,  +  \jtT  —  sin2  ^er  cos  9,  +  \jtT  —  sin2  0,J 

The  function  W(p,  q)  is  the  spatial-frequency  power  spectrum  of  the  height  of  the 
rough  surface: 

W(p,q)  =  4  I"  I"  e-ipx-'qyR(x,y)dxdy  (166) 

If  J- oo  7-00 


where  R(x,y)  =<  £(x',  t/')£(x'  —  x,  y'  —  y)  >  is  the  autocorrelation  function  of 
the  rough-surface  height  y').  The  brackets  <>  denote  the  expected  value; 
and  it  is  assumed  that  the  height  £  is  a  stationary,  zero-mean  random  process  in 
two  dimensions. 

We  shall  examine  the  special  case  in  which 


•  the  incident  wave  is  randomly  polarized; 


79 


•  the  angle  of  incidence  9{  =  0;  and 

•  the  autocorrelation  R(x,y)  is  isotropic. 

When  the  incident  wave  is  randomly  polarized,  S^nc  =  S£c  =  iSinc.  Thus 


where 


Furthermore,  when  9t  =  0  we  find 


AS* 

AA  C 

(167) 

~  4 *W<rv  ,ae 

(168) 

°h 

2^°'hh 

(169) 

°v 

2  (&vh  +  ®vv  ) 

(170) 

find 

Ohh 

=  -ah  cos  (f>3 

(171) 

Ohv 

=  ahs\n<t>s 

(172) 

Ovh 

=  —ct*,  sin  <f>. 

(173) 

ovv 

=  —av  cos  <j>. 

(174) 

where 


Oh  = 


av  — 


s/^T  ~  1 _ 

cos  9,  +  \JeT  —  sin2  0, 
(y/^  ~  l)  \/er  ~  sin2  9, 


eT  cos  6,  +  \Ji eT  —  sin2  9, 

Finally,  when  the  autocorrelation  function  can  be  written  in  the  form 

R{x,y)  =  Wf(Jx*  +  y*/l) 


(175) 


(176) 


(177) 
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where  /  is  a  characteristic  correlation  distance  and  h 2  is  the  variance  of  the  rough- 
surface  height,  we  find  that 

2/2/i2  r°° 

W(-)  =  -  /  f(u)J0(klus\n9s)udu  ( 1 7S) 

n  Jo 

=  l2h?F{kl  sin#,) 

where  J0(-)  denotes  the  Bessel  function  of  order  zero.  Thus 

ah  =  2nk4Fh2  cos2  0sF(kl  sin  03)  |  ah  |2  (179) 

<t„  =  2nk4l2h2  cos2  0,F(kl  sin  9,)  j  av  |2  (ISO) 

We  next  consider  the  power  received  by  an  antenna  situated  over  the  rough  surface. 


4.2  POWER  COLLECTED  BY  AN  ANTENNA  OVER  A  ROUGH 
SURFACE. 


Let  the  antenna  be  located  at  height  h  over  the  rough  surface  and  let  its  axial 
direction  be  defined  by  angles  0b ,  <f>B,  as  shown  in  Figure  72.  The  effective 
aperture  of  the  antenna  is  denoted  Ae(9',  <f>')  where  O'  and  <fr'  are  defined  with 
respect  to  the  antenna  axis. 

The  power  density  scattered  in  the  direction  of  the  antenna  by  a  small  area  A/1 
in  the  z  —  0  plane  is 

=  AnR2  Vh.vSinc  (181) 

and  the  power  A P  received  from  this  small  area  is 


^ph,v  _  ‘S'inc^e  (9/ >  }&h,v  A/1  COS  0S 

in  COS  9,  R2 


(182) 


Now  A/1  cos  9,/R2  is  just  AQ,  the  solid  angle  subtended  by  the  area  A  A.  Thus 
the  total  power  received  by  the  antenna  is  simply 


ph,v  _  Sine  j  Ae{ff  ^ 

4  7r  7n0  cos  9, 


(183) 
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where  denotes  the  total  solid  angle  (of  2tt  steradians)  subtended  by  the  rough 
surface. 

We  consider  the  special  case  where  the  axis  of  the  antenna  is  parallel  to  the  x-axis, 
so  that 

cos  0S  =  sin  &  cos  <f>'  (184) 

sin  9,  sin  4>a  =  sin  O'  sin  <j>'  (185) 


x 

Figure  72.  Antenna  over  a  rough  surface:  geometry. 

Furthermore,  dfi  =  sin  O'dO'dfi  in  the  antenna  coordinates  with  —  tt/ 2  <<?>'< 
7r/2,  0  <  O'  <  x.  Thus 

Ph'v  =  ~  f2  sec  W  ^  A'i?,  <f>')eUO,)dO'  (186) 

47T  J-ir/2  Jo 

Substituting  eqs.  (179),  (180),  (184),  and  (185)  into  eq.  (186),  we  obtain 

Ph-V  =  -Sinck*!2!?  r'2  cos  4>'d<f>'  [*  sin2  O'dO'  •  (187) 

2  J—ir/2  JO 
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where 


F  -  sin2  O'  cos2  <f>^J  Ae(0',<f>')  |  eth,v(Or ,  <j>) 

s/F  -  1 

ah  =  - p — = - ^  ■■ 

sin  O'  cos'  +y  er  —  1  +  sin2  O'  cos2  <f>' 

i^y/F  —  l)  \fe7~-^\^vFo^cos^' 

eT  sin  O'  cos  <j>'  -f  \Jer  —  1  +  sin2  O'  cos2  <f>' 

F{0  =  ~  [  f(u)J0(ut)udu 

7T  JO 


<y„  = 


(188) 


(189) 


(190) 


In  eqs.  (186)  and  (187)  above,  Ph,v  has  been  expressed  in  terms  of  an  integral  over 
the  antenna  coordinates  O',  ft.  It  may  also  be  written  in  terms  of  the  coordinates 
0 ,  and  4>3 .  Specifically,  defining 


Ae(0s,<f>a)  =  Ae(0',<f>') 


we  obtain 


,  v  r r*  * 

Ph’v  =  tan  03ah>v(0,)d0,  Ae(0s,  <j>s)d<t>a 

47 r  Jo  J-t 

C.  _  rir/2 

=  - ~k4l2h 2  \  Qh<v{0,)  \2  F(kl  sm  0,)  sin  0,  cos  0,d0, 

2  Jo 


(191) 


(192) 


/  A e(0„(j),)d<f>, 

J  —  7T 

This  latter  expression  may  be  more  suitable  than  that  given  in  eq.  (187)  for  a 
given  application. 


4.3  EXAMPLE  CALCULATIONS. 

It  is  evident  that  except  under  conditions  which  allow  many  simplifying  approx¬ 
imations  to  be  made,  PK'V  must  be  evaluated  numerically.  We  shall  confine  our 
attention  to  a  specific  antenna  configuration,  the  rectangular  aperture  with  uni¬ 
form  illumination.  For  such  an  antenna, 

Ae(0',<f>')  =  ab  sine2  sin#'  cos</A  sine2  ( —  sin  0'  sin  <f>'\  (193) 
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Ae(0a,  <j>3)  =  ab  sine2  cos  03j  sine2  sin  6,  sin  4>,  j  (194) 

where  a  and  b  are  the  aperture  dimensions  and  the  sides  of  length  a  are  taken  to 
be  parallel  to  the  z-axis.  Furthermore,  sinc(x)  =  (sinx)/x. 

We  shall  also  consider  that  F(£)  is  of  the  form  appropriate  for  a  Gaussian  corre¬ 
lation  function  (cf.  eq.  (177)) 

R(x,y)  (195) 

Thus 

F(  0  =  -e'*2/4  (196) 

7T 

Substituting  eqs.  (194)  and  (196)  into  eq.  (192),  we  obtain 


ph,v  _  F-,  „cab 

2tt 


sin  03  cos  0 


k4l2h2J*'2  |  Qh,v(0s )  P  e-(^sin<?')2sinc2  cos0^  •  (197) 

/*y2  '( kb  \ 

sine2  I  —  sin  03  sin  (f>3  d<f>3 

■*/ 2  \  2  / 


wherein  we  have  altered  the  limits  of  the  integral  over  <j>3  to  include  only  the 
“forward-looking”  portion  of  the  antenna  pattern. 

The  integral  over  <f>,  can  be  expressed  in  terms  of  a  function  G(£)  defined  as 


follows: 

£ 

J-n 

II 

^  1  to 

OO 

-  E 

k=0 

with 

<7(0- 

Thus 

Ph  v  =  S'm?bk4l2h*  [ 
2  Jo 

•  2  f  lea  . 

sine  —  cos  B. 
\  2 

(198) 


^  (-)*£2* 

,fe(2fc +  !)*!(*+!)! 

-*  7  (t  °°) 


(199) 

(200) 
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A  plot  of  the  function  G{£)  vs.  £  is  shown  in  Figure  73. 

We  have  performed  sample  calculations  of  the  expression  for  Ph  given  in  eq.  (200) 
for  the  case  in  which  j  er  |>>  1.  In  this  case,  =  1;  and  we  have  evaluated  the 
normalized  quantity 


Ph  = 


Ph 

abSin 


(201) 


^  r/2 


=  ~k4l2h2  f 

2  Jo 


,-(*sin  6,f 


sine 


/  ka 

u 


cos  6, 


G 


fkb  .  n\  ,  A 
2  sin  0S J  sin  0 3 


cos O.dO, 


which  is  the  ratio  of  the  power  received  from  rough-surface  scattering  to  the  total 
power  which  is  available  to  the  antenna. 

It  is  evident  from  eq.  (201)  that  the  parameters  on  which  p/,  depends  are  the 
electrical  dimensions  of  the  antenna,  ka  and  kb ;  the  electrical  correlation  length 
kl ;  and  the  electrical  rough-surface  height  variance  k2h2.  The  factor  k2h 2  can  be 
considered  as  the  limiting  form  for  small  height  variances  of  the  factor 

1  —  exp  (^—k2h2^j 

whose  maximum  value  is  unity.  In  order  to  bound  the  collected  power  we  shall 
assume  this  factor  to  be  unity  in  the  computational  results.  Thus 

Ph  <ph  =  ^ k 212  jf  e~(Wu/2)2sinc2  y/\~-  u2J  G  j  udu  (202) 


It  is  easy  to  show  that  ph  becomes,  in  the  limit  of  long  correlation  lengths  /, 


..  .  .  2ka 

hm  ph  =  sine  — 

I— *oo  Z 


(203) 
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£ 

Figure  73.  G(f)  vs.  f.  The  function  differs  significantly  from  1  /if  only  for 
0  <  f  <  3. 


This  limit  is  obviously  related  to  the  effective  area  of  the  aperture  antenna  at 
grazing  incidence  (scattering  angle  9,  =  0).  For  simplicity  in  the  numerical  results 
to  follow  we  shall  take  ka  =  2 ns  (n  >  1).  Thus  the  limit  of  f>h  for  large  /  will  be 
zero  in  every  case  which  we  consider.  The  maximum  gain  Go  of  the  aperture  is 


1  b  2  4n27r6 

G0  =  --(ka)"  =  - 

7 r  a  a 


(204) 


We  shall  choose  ratios  bja  such  that  (204)  is  satisfied  for  a  given  value  of  Go- 
Thus  ph  can  be  written 


(  mrl  \  2 

[  ( nirl  \  21 

-  /  exp 

1 

1 

\  a  J  Jo 

L  v  «  /  j 
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wherein  Go,n  and  l/a  are  parameters.  We  have  evaluated  as  a  function  of  l/a 
for  the  cases  listed  below: 


Go 

=  10 

dB, 

n 

=  1 

(b/a  = 

=  0.80) 

Go 

=  15 

dB, 

n 

=  1 

(b/a  = 

=  2.52) 

n 

=  2 

(b/a  = 

=  0.63) 

Go 

=  20 

dB, 

n 

=  2 

(b/a  = 

=  1.99) 

n 

=  3 

(6/a  = 

=  0.88) 

n 

=  4 

(b/a  -- 

=  0.50) 

-10 


Figure  74.  Normalized  received  power  ph.  vs.  normalized  correlation  length  l/a. 
G0  is  given  in  dB. 
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The  results  are  shown  in  Figure  74.  It  is  evident  from  the  curves  shown  there  that 


1.  The  peak  value  of  ph  as  a  function  of  l/a  decreases  as  the  antenna  gain  is 
increased  and  as  n  is  increased. 

2.  The  value  of  l/a  at  which  the  peak  occurs  decreases  as  the  gain  is  increased 
and  as  n  is  increased,  and  is  less  than  0.5  in  all  curves  considered. 

3.  The  peak  value  of  ph  in  dB  is  less  than  the  negative  of  the  antenna  gain  in 
dB. 

4.  The  peak  value  of  p^  is  not  a  strong  function  of  n. 

Results  3  and  4  indicate  that  at  least  for  the  configuration  considered  here,  the 
maximum  power  collected  from  rough  surface  scattering  is  easily  estimated  and 
that  it  is  not  a  strong  function  of  the  antenna  pattern. 

4.4  CONCLUDING  REMARKS. 

We  have  considered  the  scattered  power  collected  by  an  antenna  over  a  rough 
surface  and  have  calculated  some  representative  results  for  a  rectangular  aperture 
antenna  whose  axis  is  directed  parallel  to  the  surface.  It  has  been  shown  that 
the  collected  power  decreases  as  the  antenna  gain  increases  and  that  this  power 
is  maximized  when  the  coherence  length  of  the  scattering  surface  is  less  than  the 
antenna’s  vertical  dimension.  The  peak  value  of  the  collected  power,  normalized 
to  the  power  available  to  the  antenna,  is  roughly  equal  to  the  negative  of  the 
antenna  gain  in  dB  and  does  not  depend  strongly  on  the  details  of  the  antenna 
pattern. 

Our  treatment  of  this  problem  has  been  based  upon  an  analytical  model  derived 
from  a  first-order  perturbation  theory.  This  theory  provides  the  scattering  cross- 
sections  <7hhi ...  in  terms  of  the  properties  of  the  scattering  medium  and  its  rough 
surface.  If  experimental  data  for  these  quantities  are  used,  the  framework  of  the 
analysis  given  herein  remains  unchanged.  The  result  given  in  (183)  would  simply 
employ  experimentally  derived  expressions  for  <7h  and  crv  rather  than  those  used 
here. 
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SECTION  5 

APERTURE  COUPLING  TO  A  WIRE  IN  A  CYLINDRICAL 

CAVITY 


Summary:  We  examine  the  electromagnetic  coupling  to  a  thin  wire  on  the  axis  of 
a  circular  cylindrical  cavity  through  a  thin  azimuthal  slot  aperture  on  the  cavity 
side  wall.  The  wire  is  shorted  to  a  cavity  end  wall  at  one  end  and  it  becomes  the 
center  conductor  of  a  coaxial  cable  at  the  other  end.  The  coupling  is  described  in 
terms  of  the  Thevenin  equivalent  circuit  driving  the  feed  point  of  the  coaxial  cable, 
so  that  an  arbitrary  load  may  be  considered.  The  source  is  specified  in  terms  of  the 
voltage  across  the  slot  aperture  in  the  outer  wall,  so  that  the  exterior-to-interior 
coupling  problem  may  be  approximately  decoupled.  Representative  numerical 
data  are  presented  to  illustrate  the  results. 

In  this  section  we  examine  a  relatively  simple  and  analytically  tractable  approxi¬ 
mate  treatment  of  electromagnetic  coupling  through  an  aperture  in  a  cavity  wall 
to  a  wire  inside  the  cavity.  Our  purpose  is  to  develop  the  transfer  function  which 
relates  the  incident  electromagnetic  field  on  the  exterior  of  the  cavity  to  the  cur¬ 
rent  induced  in  an  internal  load. 

The  geometry  of  the  problem  is  shown  in  Figure  75.  The  cavity  itself  is  a  right 
circular  cylinder  of  radius  b  and  length  /.  A  wire  radius  of  a  is  coaxial  with  the 
cavity  and  shorted  to  the  cavity  wall  at  r  =  0.  It  leaves  the  cavity  at  z  =  /, 
forming  the  center  conductor  of  a  coaxial  line  whose  outer  radius  is  c.  We  assume 
that  a  <<  b  and  c  <<  6.  The  cavity  is  excited  by  a  thin  circumferentially-oriented 
slot  of  width  w  and  angular  extent  A<£,  located  at  z  =  Zq  in  the  cavity  side  wall. 

In  our  analysis  of  the  exterior  problem  the  slot  is  treated  as  if  it  were  located  in 
an  infinitely  extended  conducting  plane.  Use  is  made  of  results  in  Section  2  to 
estimate  the  voltage  across  the  slot.  This  slot  voltage  is  then  used  to  drive  the 
interior  problem.  Thus  the  transfer  function  from  the  incident  field  to  the  internal 
wire  current  is  the  product  of 

•  the  transfer  function  relating  the  incident  field  to  the  slot  voltage,  derived 
in  Section  2;  and 

•  the  transfer  function  relating  the  slot  voltage  to  the  internal  wire  current, 
derived  herein. 
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The  approximation  employed  herein  obviously  does  not  take  into  account  the 
effect  of  the  cavity  itself  on  the  fields  in  the  slot.  This  approximation  will  cause 
the  cavity  Q  to  appear  higher  than  it  actually  is,  because  reradiation  from  the 
slot  is  not  accounted  for. 


Z=0  Zo  £ 


Figure  75.  Geometry  of  the  problem.  The  slot  width  is  w  and  the  outer  radius 
of  the  coax  leaving  the  cavity  is  c. 


5.1  FORMULATION. 


The  electromagnetic  field  which  excites  axial  currents  on  the  center  conductor  in 
the  cavity  is  the  axially  symmetric  TMZ  portion  of  the  total  field,  expressed  in 
terms  of  a  scalar  function  0  by  the  relations 

B  =  J-  ft. 

p  ju>e0  dpdz 


Ez 


1 

jut  o 


0 


H * 


d\l> 

dp 
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The  function  rp  is  a  solution  of  the  scalar  Helmholtz  equation  in  circular-cylindrical 
coordinates 


dz 2 


+  k?xp  =  0 


(207) 


wherein  k 2  =  u>2fi0e o  with  (i0  and  e0  denoting  the  permeability  and  permittivity  of 
free  space  respectively.  The  time  dependence  exp  {jut)  is  assumed  and  suppressed. 


Our  approach  to  solving  the  problem  of  electromagnetic  coupling  from  a  slot 
aperture  on  the  outer  side  wall  of  the  cavity  to  a  transmission  line  which  leaves 
the  center  of  an  end  wall  is  as  follows:  first,  the  current  70  on  the  center  conductor 
due  to  the  excitation  of  the  outer  slot  is  determined  as  a  function  of  the  average 
slot  voltage  V0  and  written  in  the  form 


i0  =  v0v0  (Vi  =  0)  (208) 

In  the  calculation  of  70  it  is  assumed  that  the  coupling  aperture  to  the  transmission 
line  is  closed.  Next,  the  current  7t  due  to  voltage  V)  across  the  transmission-line 
aperture  is  determined,  assuming  that  the  outer-wall  slot  is  closed,  so  that 

h  =  Vi  Vi  (Vo  =  0)  (209) 

The  total  current,  70-f  7i,  is  then  related  to  the  voltage  Vi  by  the  input  admittance 
to  the  transmission  line  Vi: 

70  +  7,  =  YiVi  (210) 

The  relation  (210)  is  used  to  solve  for  Vi  in  terms  of  V0.  Thus 

V1  =  Y0%(Yi-Yl)~l  (211) 

and  the  solution  is  complete.  Knowledge  of  the  external  excitation  (Vo)  and  the 
line  input  admittance  (Vj)  yields  the  line  excitation  (Vi).  In  the  following,  we 
carry  out  each  of  the  steps  leading  to  the  result  (211). 

5.1.1  Coupling  through  the  Outer  Slot. 

The  electric  field  Ez  at  the  outer  slot  at  p  —  6,  |  <f>  |<  A<j!>/2,  z  =  zq  is  presumed 
known: 

Ez,  =  V0(| i>)6{z  -  Zq)  ( p  =  b)  (212) 
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We  assume  that  E#  =  0  over  the  slot.  The  function  rp  is  conveniently  represented 
in  the  form 

oo  /  Y17T  z\  1 

xf>  =  4> 0  =  jut o  Y  COS  (  — —  )  —  Co(Anp,  A na)  (213) 

where  the  coefficients  An  are  to  be  determined  and 

C0(x,  y)  =  J0(x)N0(y)  -  J0(y)N0(x)  (214) 

A 2n  =  k2  -  (nir/l)2  (215) 

The  form  for  ip0  given  by  (213)  and  (214)  above  ensures  that  Ep  vanishes  at  z  =  0 
and  z  =  /  and  that  Ez  vanishes  at  p  =  a\  Ej,  is  identically  zero  throughout  the 
cavity. 

Matching  the  boundary  condition  (212)  in  the  least-squares  sense,  we  readily 
obtain  for  An 


A0  = 

An  = 


Vo 


lCo(kb,  ka) 

2V0  cos(mrz0/l) 

A na) 


(216) 


(n  >  1) 


where 


_  1  r^4>/2 

V0=—  V0(<f))d<f> 

C'K  J-&<1>/2 


The  current  on  the  center  wire  at  the  endpoint  z  =  l  is  given  by 

7°  =  ilu,tolVa{wc^HM) 

,  2  y'  (~)n  cos(n7rz0//) ) 

A^/2C0(An6,  Ana)  j 

=  Y0% 

which  completes  the  solution  of  the  first  part  of  the  problem. 


(217) 


(218) 
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5.1.2  Current  Due  to  End-Slot  Excitation. 
A  suitable  representation  for  the  function  ip  is 
ip  =  ip x  =  jcotoBo  In  —  cos  kz 


(219) 


+  2 juto  53  Bn  cos  (yjk2  -  hi  z)  C0(hnp,  hna ) 

n=l  '  1 

where  the  coefficients  hn  are  obtained  from  the  roots  of  the  equation 

C0{hnb,  hna)  =  0  (220) 

The  form  chosen  for  ipi  ensures  that  Ez  vanishes  at  p  =  a  and  p  =  b  and  that  Ep 
vanishes  at  2  =  0. 

The  radial  electric  field  Ep  at  z  ~  l  is  given  by 

£.(»  =  *)=  7  \  .  (a<p<c)  (221) 

p\n(c/a) 

where  c  <  b.  Vi  is  the  voltage  across  the  coaxial  line  leaving  the  cavity.  Imposing 
the  boundary  condition  (221)  leads  to  the  relation 

-  —  sin  kl  -  2  £  Bnyjk 2  -  h2n  sin  Uk7  -  h2nl)  ■  (222) 

P  n=l  '  ' 

v, 

h'n(-'o{b'np,  hnQ. )  —  J  z  z  r  {a  5;  P  —  *-) 
pln(c/a) 

wherein 

C'0(x,y)  =  Jo(x)N0(y)  -  J0(y)K(x)  (223) 

=  J0{y)Ni(x)  -  Ji{x)N0{y) 


=  -Ci(x,y) 

Integrating  (222)  with  respect  to  p  from  p  =  a  to  p  =  6  and  using  (220)  yields 

~V\ 


Bo  — 


k sin  kl\n(b/a) 


(224) 


93 


Then  multiplying  (222)  through  by  pC'0(hmp ,  hma )  and  integrating  with  respect 
to  p ,  we  obtain 


R _ -ViC0(hmc,  hma ) _ 

2^/P  -  hi  sin  (y/k*-  h2ml)  In {c/a)N„ 


(225) 


rhmb  2 

Nm  — *  /  u  [C0(ti,  ^m^)]  du 


(226) 


h2  b 2  9 

=  ^~c?(hmb,  hma)  -  4 

2  7T‘ 

Now  the  current  7t  on  the  center  conductor  at  z  =  /  is  given  by 


/j  =  4jweo/Vj 


7T  cot  kl 

2kl\n{b/a) 


(227) 


o  ^  cot  [\Jk2  -h2  /)  C0(hnc,  hna ) 
^  2t/fc2  -  hi  lNn  In  (c/a) 


=  YXVX 


We  shall  assume  that  both  c  and  a  are  small  and  make  use  of  the  fact  that  in  this 

limit  ^ 

Co(hnc,  hna)  = - In  -  (228) 

7r  a 


Thus 


Ki  =  2njuitol 


cot  kl 
kl\n(b/a) 


(229) 


^2, _ cot  [yjk2  -  h2  /) _ 

+  hx  (yjk2  -h2  l)  \(*hnb/2)2  C2(hnb,  hna)  -  1 
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5.1.3  Exterior  to  Coaxial  Line  Coupling. 

The  total  current  on  the  center  conductor  at  the  feed  point  of  the  coaxial  cable 
is  just  To  -f-  I\.  The  two  current  components  are  now  related  by  imposing  the 
condition  that  the  coaxial  cable  leaving  the  cavity  presents  an  impedance  Z,  at 
the  point  z  —  l.  Thus,  with  Fj  =  1/Z*,  we  have 


/o  +  A  =  Y,Vy  (230) 

Using  (218)  and  (227)  and  solving  for  Vj,  we  obtain 

Vx  =  Y0(Yt-Y,r'V0  (231) 

which  is  the  desired  relation  between  the  “external”  excitation  represented  by  V0 
and  the  “internal”  excitation  V\. 

We  may  also  construct  an  equivalent  circuit  representing  the  source  at  the  feed 
point  of  the  coaxial  cable.  The  open-circuit  voltage  at  the  feed  point  is 


Voc  =  -YoVo/Y,  (232) 

and  the  short-circuit  current  is 

Ac  —  Iq  —  YoVo  (233) 

Thus  the  equivalent  circuit  takes  either  of  the  forms  shown  in  Figure  76.  Such 
equivalent  circuits  can  be  used  to  determine  the  voltages  and  currents  induced  in 
an  arbitrary  load. 

5.2  REPRESENTATIVE  NUMERICAL  RESULTS. 

We  present  numerical  results  for  the  transfer  functions  relating  the  current  or 
voltage  induced  in  a  50  Q,  load  to  the  average  slot  voltage  Vo  and  to  an  incident 
electric  field. 
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-wx 


(a)  Th6venin  (b)  Norton 


Figure  76.  Equivalent  circuits  for  the  feed  point  of  the  coaxial  line. 

The  current  in  a  resistive  load  Rp  connected  across  the  terminals  of  either  of  the 
equivalent  circuits  shown  in  Figure  76  is  easily  shown  to  be 


,  VoVo 

L  RLY\  -  1 

(234) 

We  define  a  transfer  admittance  Yp  as 

„  II  Vo 

T  Vo  RLYt  - 1 

(235) 

We  have  calculated  Yp  as  a  function  of  frequency  for  the 

“baseline”  configuration: 

Rl  =  50  fi 
/  =  40  cm 

z0  =  20  cm 

6  =10  cm 

a  =  0.5  mm 

(load  resistance) 
(cylinder  length) 
(aperture  position) 
(cylinder  radius) 

(wire  radius) 

as  well  as  excursions  from  this  baseline,  over  the  frequency  range  0-6  GHz.  Plots 
of  the  magnitude  of  Yp  are  shown  in  Figures  77  to  84. 
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frequency  (GHz) 


Figure  78. 


Magnitude  of  Yj  vs.  frequency  for  altered  configuration:  / 
30  cm,  b  =  10  cm,  z0  =  20  cm,  a  =  0.5  mm,  Ri  —  50  ft; 


0<  /  <  6  GHz. 


Figure  79. 


11 

frequency  (GHz) 

Magnitude  of  Yp  vs.  frequency  for  altered  configuration:  / 
50  cm,  6=10  cm,  z0  =  20  cm,  a  =  0.5  mm,  ftp  =  50  ft; 
0<  /  <  6  GHz. 


1.0 


2.0 


3.0  1.0  5.0 

frequency  (GHz) 

;  of  Yt  vs.  frequency  for  altered  configuration:  / 
=  10  cm,  z0  =  10  cm,  a  =  0.5  mm,  Ri  —  50  Q; 
GHz. 
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Figure  81.  Magnitude  of  Yj  vs.  frequency  for  altered  configuration: 

l  =  40  cm,  6  =  10  cm,  zQ  =  30  cm,  a  =  0.5  mm,  Ri  = 
50  fi;  0<  /  <  6  GHz. 
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frequency  (GHz) 


Figure  82.  Magnitude  of  Yt  vs.  frequency  for  altered  configuration: 

/  =  40  cm,  6  =  5  cm,  zQ  =  20  cm,  a  =  0.5  mm,  RL  =  50  0; 
0<  /  <  6  GHz. 


frequency  (GHz) 


Figure  84.  Magnitude  of  Yj  vs.  frequency  for  altered  configuration 
/  =  40  cm,  b  =  10  cm,  z0  =  20  cm,  a  =  1  mm,  RL  =  5i 
0<  /  <  6  GHz. 


The  highly  resonant  character  of  the  configuration  is  evident  in  each  of  these  plots. 
The  TEM  mode  resonances  are  seen  at  the  lowest  frequencies  and  are  especiallv 
clear  in  Figure  82,  where  the  cavity  radius  was  chosen  to  be  5  cm.  The  higher- 
order  mode  resonances  generally  cause  |  Vj  |  to  peak  at  higher  values  than  the 
TEM  mode  resonances  and  the  Q  of  the  higher-order  modes  is  also  higher.  The 
peak  values  of  |  Yj  |  in  every  case  are  a  few  tenths  of  a  mho.  For  a  load  impedance 
of  50  Cl,  therefore,  the  load  voltage  at  the  resonant  peaks  is  of  the  order  of  5  to 
10  times  Vo,  the  average  slot  voltage.  Away  from  the  resonances,  the  load  voltage 
is  of  the  order  of  the  average  slot  voltage  or  less. 

It  is  of  greater  practical  interest  to  consider  the  exterior-to-interior  coupling  prob¬ 
lem  in  its  entirety.  To  do  so,  we  need  to  >•••  v  the  average  slot  voltage  Vo  in  terms 
of  the  field  incident  upon  the  cylind(  from  the  outside.  We  may  use  the  results  of 
Section  2  to  show  that  the  voltage  across  a  slot  of  width  w  in  an  infinitely  extended 
conducting  plane  illuminated  by  a  normally  incident  plane  electromagnetic  wave 
whose  electric  vector  is  dired^d  across  the  slot  is  given  by 

2njEj  (  f  A: 

V (z  )  =  - —  (  sin  kl  —  2  sin  — 

ClD(kl)  V  2 


cosArz'j  (|  *'  |<  i/2) 


(236) 


wherein  l  denotes  the  slot  length,  z'  is  a  coordinate  along  the  long  dimension  of 
the  slot,  Cl  is  the  antenna  parameter 


Cl  =  2  In 


and 


1  _r  1  ,  l  ,  ^ 

.h  n.ho  ’’  +  2s 


D(ki)  (ki)2  6  nfrc 


717T 


1 


+ 


1 


kl  —  nir  —  jan  n7r  +  ;anJ 


The  damping  constants  an  are  given  by 


1 


=  [7  +  ln(2|n  1 7r )  -  Ci(2 1  n  j  tt)] 


(237) 


(238) 


(239) 


wherein  7  =  0.57721  ...  (Euler’s  constant)  and  Ci(  )  denotes  the  cosine  integral 
function  [13]. 

‘The  prime  on  the  summation  sign  indicates  that  the  n  =  0  term  is  to  be  omitted. 
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We  shall  use  the  expression  (236)  to  describe  the  voltage  across  the  slot  in  the 
cylindrical  surface,  assuming  that  the  total  angular  extent  of  the  slot  A <j>  is  less, 
say,  than  one  radian.  The  average  slot  voltage  is  easily  shown  to  be 


jEj2 

bQD(ki) 


sin  kl - :  sin 

kl 


2 


kl 

2 


(240) 


We  show  in  Figures  85  and  86  plots  of  the  magnitude  of  the  quantity 


27r6— -  2  xj 

■  vo  — 


EiP  QD(kE 


—  (sin  ki - :  sin2  — 

/)  V  kl  2  , 


(241) 


as  a  function  of  kl  (0  <  kl  <  2r)  for  ft  =  5  and  f l  —  10. 

The  voltage  Vi  across  a  load  resistance  Ri  driven  by  the  equivalent  sources  in 
Figure  76  can  be  expressed  in  terms  of  Vo  and  the  transfer  admittance  Yt  as 


VL  =  V0YtRl 


jE,i 2 

bnD{ki) 


YqRl 

yxrl-  i 


(242) 


We  have  calculated  the  transfer  function  V£,/(£,7)  as  a  function  of  frequency  for 
the  baseline  configuration  and  excursions  described  earlier  in  this  section,  under 
the  assumption  that  the  slot  length  /  =  5  cm.  The  fundamental  slot  resonance 
occurs  at  3  GHz.  Plots  of  the  magnitude  of  this  transfer  function  as  a  function  of 
frequency  are  shown  in  Figures  87  -  94  for  fi  =  5  and  95  -  102  for  fl  =  10. 

The  broadly  resonant  character  of  the  average  slot  voltage  as  a  function  of  kl  is 
evident  in  the  curves  shown  in  Figures  85  and  86.  The  peak  value  of  the  quantity 
plotted  there  is  almost  unaffected  by  the  slot  width  w,  but  the  sharpness  of  the 
resonance  is  increased  slightly  as  w  is  decreased.  The  overall  effect  of  the  slot 
response  on  the  normalized  load  voltage  Vl/(E,7)  is  simply  to  “envelope”  the 
many  cavity  resonances,  as  shown  in  Figures  87  -  102.  The  peak  values  of  the 
magnitude  of  the  load  voltage  are  in  every  case  less  than  Ej  and  are  typically  only 
a  few  tenths  of  Eil.  It  is  easy  to  show  that  the  effective  area  of  the  slot  penetration 

A 

is  thus  of  order  r  at  a  resonance  of  the  interior  cavity  and  considerably  less  away 
from  resonance. 
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FREQUENCY  (GHZ) 


Figure  87.  Magnitude  of  Vl,/(EJ)  vs.  frequency  for  baseline  configuration: 

(  =  40  cm,  6  =  10  cm,  z0  =  20  cm,  a  =  0.5  mm,  Ri  =  50  ft, 
l  =  5  cm;  ft  =  5. 


FREQUENCY  (GHZ) 

Figure  93.  Magnitude  of  Vl/(E, l)  vs.  frequency  for  baseline  configuration, 
except  6  =  15  cm;  ft  =  5. 
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0. 


Figure  95.  Magnitude  of  Vi/(Eil)  vs.  frequency  for  baseline  configuration: 

(  =  40  cm,  b  =  10  cm,  zQ  =  20  cm,  a  =  0.5  mm,  RL  =  50  Q, 
l  =  5  cm;  ft  =  10. 
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5.3  CONCLUDING  REMARKS. 

We  have  considered  a  simple  problem  of  microwave  coupling  to  a  load  in  the 
interior  of  a  shielded  region,  the  penetration  occurring  via  a  thin  slot  aperture 
oriented  and  excited  so  as  to  yield  the  maximum  coupling  to  the  load.  The 
treatment  of  the  problem  was  approximate  in  the  sense  that  the  effects  of  the 
cavity  on  the  average  slot  voltage  were  not  taken  into  account.  This  approximation 
causes  the  cavity  resonances  to  appear  sharper  than  they  actually  are,  because 
the  reradiation  of  the  cavity  fields  out  of  the  slot  is  not  accounted  for. 

For  the  configuration  which  was  explored  numerically,  the  cavity  resonances  were 
much  sharper  than  that  of  the  slot  aperture,  so  that  the  frequency-domain  mag¬ 
nitude  of  the  load  voltage  (normalized  to  an  incident  voltage  Ej)  was  dominated 
in  its  structure  by  the  sharp  cavity  resonances.  The  effective  area  of  the  slot  aper¬ 
ture  was  of  the  order  of  the  square  of  the  slot  length  at  a  resonance  of  the  cavity, 
decreasing  to  zero  at  anti-resonances.  The  effective  area  for  a  pulsed  signal  will 
depend  upon  the  center  frequency  and  bandwidth  of  the  pulse. 

Several  extensions  and  variants  of  the  present  problem  would  be  of  interest  for 
microwave  applications.  One  such  variant  would  be  the  problem  of  coupling  to  a 
cavity  homogeneously  filled  with  an  absorptive  material  representing,  say,  circuit 
boards  and  other  electronic  components.  One  would  compute  the  total  energy 
absorbed  by  the  lossy  material  under  pulse  excitation  of  the  cavity.  It  would  also 
be  a  worthwhile  exercise  to  refine  the  treatment  of  the  problem  treated  herein  to 
include  the  effects  of  the  cavity  fields  upon  the  slot  voltage. 

Microwave  coupling  applications  of  interest  at  the  present  time  will  tend  to  involve 
cavities  which  are  reasonably  large  in  comparison  to  the  wavelength,  and  coupling 
apertures  which  are  comparable  to  the  wavelength  in  size.  Thus  the  number  of 
cavity  resonances  “spanned”  by  the  aperture  response  will,  as  in  the  case  treated 
in  this  section,  be  large. 

It  is  an  interesting  exercise  to  estimate  the  total  power  or  energy  coupled  to  a  load 
inside  a  general  slot-excited  cavity.  It  is  easily  shown  that  the  number  of  cavity 
modes  having  resonance  frequencies  between  /  and  /  +  A/  is  approximately 

AiV  =  -/2A/  (243) 

c5 

where  V  is  the  cavity  volume  and  c  denotes  the  speed  of  light.  A  signal  of  center 
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frequency  f0  and  bandwidth  A f0  will  thus  excite  approximately 

AJV„  = 


(244) 


modes.  Denote  the  quality  factor  of  these  (loaded)  resonant  modes  by  Q 2.  Thus 
the  width  A/r  of  a  given  resonance  will  be  approximately 


A  fr  =  £ 

Q 

and  the  power  absorbed  in  a  single  resonance  will  be 

p  ~  c.  a  — 

1  r  —  Jinc^e  « 


(245) 


(246) 


where  Smc  denotes  the  incident  power  density  spectrum  (watts  per  square  meter 
per  Hz)  and  Ae  is  an  effective  area.  The  total  power  absorbed  Pa  is  thus 


Pa  =  PtAN0 


(247) 


—  ^inc^/o^i 


SttV/o3 


—  •Winc'd 


8  ttT/3 


where  5jnc  =  SmcAfo  is  the  incident  power  density  (watts  per  square  meter). 
Now  if  the  center  frequency  fo  is  that  of  the  (broad)  slot  resonance,  we  have 


(248) 


Furthermore,  if  Ae  at  a  cavity  resonance  can  be  written  as 


where  a  is  a  constant,  we  have 


A,  =  or 


P  Z-S- 

‘  Qi  " 


(249) 


(250) 


2It  is,  of  course,  an  approximation  to  ascribe  the  same  value  for  Q  to  each  of  the  (many) 
modes  under  consideration. 


Thus  the  effective  area  for  power  absorbed  from  the  incident  signal  is 


anV 

A*  =  —f  (251) 

Vt 

where  the  center  frequency  is  close  to  that  of  the  slot  resonance  and  when  the 
bandwidth  of  the  incident  signal  includes  several  cavity  resonances.  For  the  ge¬ 
ometry  considered  in  this  chapter,  and  estimating  Q  =  50  and  a  =  l/ff,  we  find 
Aeff  =  16tt  cm2. 
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SECTION  6 

WAVEGUIDE  FIELD  COUPLING  TO  A  FILAMENTARY  LOAD 


Summary:  A  model  is  constructed  to  represent  the  Thevenin  equivalent  circuit 
which  drives  a  thin  coaxial  probe  in  a  rectangular  waveguide.  The  model  is  in¬ 
tended  to  represent  the  drive  on  a  detector/mixer  diode  in  a  microwave  receiver. 
Representative  numerical  results  are  presented  which  display  the  coupling  to  a 
linear  resistive  load  as  a  function  of  the  various  problem  parameters,  over  the  fre¬ 
quency  band  in  which  only  the  dominant  mode  can  propagate  in  the  waveguide. 
It  is  found  that  even  when  the  probe  reactance  is  uncompensated,  the  absorbed 
power  can  still  reach  within  3  dB  of  the  incident  power  for  the  example  cases 
considered. 

An  important  problem  in  microwave  system  vulnerability  analysis  is  that  of  de¬ 
termining  the  coupling  of  energy  which  has  entered  the  system  to  a  potentially 
vulnerable  detector  or  mixer  -  typically  a  very  sensitive  diode.  The  diode  is 
typically  coupled  to  the  electromagnetic  field  in  the  waveguide  by  a  thin  coaxial 
probe  which  extends  into  the  guide.  The  detected  signal  is  then  observed  across 
the  diode.  We  consider  a  model  for  the  equivalent  circuit  which  drives  the  diode. 
We  assume  that  the  waveguide  is  terminated  beyond  the  coaxial  probe  location 
by  a  short  circuit.  This  configuration  is  typical;  the  position  of  the  short  is  ad¬ 
justed  to  maximize  the  signal  at  the  probe  position,  at  the  center  of  the  operating 
frequency  band. 

In  the  remainder  of  this  section  we  set  up  an  equivalent  circuit  at  the  probe  ter¬ 
minals  and  evaluate  the  elements  in  this  circuit.  We  present  numerical  results 
which  show  the  power  coupled  to  a  linear  load  representing  the  small-signal  diode 
impedance.  In  particular,  we  explore  the  frequency  dependence  of  the  absorbed 
power  over  the  operating  frequency  band.  Implications  for  vulnerability  assess¬ 
ment  are  discussed. 

6.1  EQUIVALENT  CIRCUIT  FOR  A  WAVEGUIDE  PROBE. 

The  geometry  of  the  waveguide/probe  configuration  is  shown  in  Figure  103.  The 
waveguide  is  rectangular  in  cross-section,  of  dimensions  axb  (a  >  b),  and  extends 
from  z  =  — oo  to  z  =  /,  where  it  is  terminated  in  a  short  circuit.  The  probe  is 
the  center  conductor  of  a  coaxial  connector  on  the  bottom  ( y  =  0)  wall  of  the 
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waveguide;  its  radius  is  denoted  r  and  it  is  centered  at  ar  =  x0,  z  ■=  0.  It  is 
assumed  that  r  «  a  and  that  the  frequency  of  operation  is  such  that  only  the 
dominant  TEi0  mode  can  propagate  in  the  waveguide. 

A  transmission-line  equivalent  circuit  for  this  configuration  is  shown  in  Figure 
104  [14].  The  purpose  of  the  analysis  of  the  next  subsection  is  to  evaluate  the 
transformer  ratio  N  and  the  reactance  Xp.  Before  proceeding,  however,  we  reduce 
the  transmission-line  model  of  Figure  104  to  a  Thevenin  equivalent  circuit  at  the 
terminals  of  the  coaxial  connector. 

We  assume  that  a  TE1()  mode  field  is  incident  on  the  probe  and  the  waveguide 
termination  from  z  —  —  oo.  The  incident  electric  field  across  the  center  of  the 
guide  will  be  denoted  E0,  so  that  the  incident  signal  voltage  Vo  =  Eob  (recall  that 
b  is  the  height  of  the  waveguide).  The  characteristic  impedance  of  the  TE]0  mode 
will  be  denoted  Z0.  Transforming  the  source  and  termination  to  the  probe  location 
yields  the  equivalent  circuit  shown  in  Figure  105.  The  transformed  termination 
impedance  Z'T  is 

Z't  —  j Zo  tan  /?io/  =:  (252) 

where  /3l0  denotes  the  propagation  constant  of  the  dominant  mode. 


side  view 


end  view 


Figure  103.  Geometry  of  the  waveguide  and  probe  configuration. 
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Once  the  parameters  N  and  Xp  have  been  determined,  the  circuit  model  can  be 
exercised  with  an  arbitrary  load.  A  diode  load  will,  of  course,  be  nonlinear;  in 
that  situation  a  circuit-analysis  code  would  be  used  for  accurate  computation  of 
the  diode  voltage  and  the  power  dissipated  in  the  diode.  However,  it  is  instructive 
to  consider  a  linearized  diode  model.  We  describe  the  diode  as  a  pure  resistance 
Rd.  Then  the  power  absorbed  in  the  diode,  normalized  to  the  incident  power,  is 
g.vcn  by  the  expression 


aba 


=  4N2RdZ0 


jXi 


JXT(Rd  +  jXp)  +  Z0(Rd  +  jXp  +  jN2XT) 


(253) 


where  the  incident  power  is  expressed  as  1 

Pin,  =  X.  |  E0  I2  (254) 

We  shall  numerically  explore  this  normalized  power  later.  We  turn  now  to  the 
evaluation  of  the  transformer  ratio  N  and  the  reactance  Xp. 


6.2  PROBE  INPUT  IMPEDANCE. 

Since  the  parameters  N  and  Xp  are  independent  of  the  details  of  the  waveguide 
termination  (a  shorting  plate  at  z  =  l)  we  address  the  computation  of  the  probe  in¬ 
put  impedance  by  considering  the  waveguide  to  be  matched  at  z  =  ±oo.  Further¬ 
more,  although  we  can  choose  a  variety  of  probe  configurations  in  the  waveguide, 
we  consider  only  the  case  where  the  probe  runs  all  the  way  across  the  waveguide 
and  is  connected  to  the  wall  at  y  =  b.  The  modeling  of  other  probe  configurations 
can  be  carried  out  using  methods  employed  in  this  section. 

We  assume  that  the  current  on  the  probe  is  of  the  approximate  functional  form 

I(y )  =  70sec  feb cos  k(y  —  6)  (255) 

The  electromagnetic  field  in  the  waveguide  is  TMy  and  derivable  from  a  scalar 
function  ip  via  the  relations 

H  =  V  x  t pay  (256) 

'This  expression  for  the  incident  power  results  from  the  choice  Vo  =  Eob  in  the  Thevenin 
equivalent  circuit. 
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E  =  — V  xVxK  (257) 

]ueQ 

where  x /»  satisfies  the  scalar  Helmholtz  equation 

V2i/>  +  k2xP  =  0  (258) 

The  propagation  constant  &  =  u>/c,  where  denotes  the  frequency  (the  frequency 
dependence  exp(jcot)  is  assumed  and  suppressed)  and  c  is  the  speed  of  light. 

We  first  model  the  surface  density  Js  in  the  plane  2  =  0  as  resulting  from  a 
zero-thickness  filament.  Thus 


Js  =  70sec  kb  cos  k(y  —  b)8(x  —  xQ)ay 
Writing  the  function  xp  as 


m=l  n= 0 


.  rmcx  niry  , 

sin - cos  exp {-jpmn  2 

a  b 


(259) 


(260) 


wherein 

Pmn  =  k2  ~  (mir/a)2  -  (nx/b)2  (261) 

and  employing  (256)  and  (259),  we  write  the  boundary  condition  at  2  =  0  as 

TTiTcx  mry 

2 j  V)  V)  PmnO-mn  sin - cos  — —  =  I0  sec  kb  cos  k(y  —  b)6(x  —  x0)  (262) 

m=l  n=0  a  ® 

The  coefficients  amn  are  then  easily  shown  to  be 

V„  .  twtxq  fb  mry 

amn  =  ~a - 7 - 77  sin -  /  COS  k(y  -  6)  cos  — r — rfj/  (263) 

jpmnabcos  kb  a  Jo  b 

/c/oentanfc6  sin(m7rxo/a) 
j0mnab  [fc2  -  {mr/b)2) 

where  Neumann’s  number  <=.„  =  1  if  n  =  0  and  en  =  2  if  n  >  1. 

We  evaluate  the  complex  power  supplied  by  the  probe  current  by  performing  the 
integral 

P  =  —  -  f  Ey-^—Io  sec  £6  cos  fc(y  —  b)dSp  (264) 

2  J Sp  2?rr 
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where  Sp  denotes  the  probe  surface.  Since  we  have  assumed  that  the  probe  ra¬ 
dius  is  small,  Ey  will  be  essentially  uniform  around  the  probe  so  that  P  is  given 
approximately  by 

/*  fb 

P  =  --~  Ey(x0,  y,r)  sec  kbcos  k(y  —  b)dy  (265) 

2  Jo 

Now  since  P  can  also  be  written  in  terms  of  the  input  impedance  Z,n  seen  at  the 
probe  terminals  from  the  coaxial  connector  side,  we  have 

P  =  \zin  |  70  |2  (266) 

Thus  using  (265)  and  (266),  we  find 


Zin  =  j^Hob 


/tan  kb\2  ^2,  ~  en  sin2 (rmr x0/ a) e~^mnT 
V  kb  )  ™=1»=0  jpmna\\  -  {rnijkby) 


(267) 


On  the  waveguide  side  of  the  transformer  in  Figure  104,  the  input  impedance  is, 
by  inspection,  Z0/2.  2  The  real  part  of  Z,„  in  (267)  must  transform  to  Z0/2;  thus 

N2  =  ~^(Zin)  (268) 


sin 


2  7T Xq 


where  we  have  used  the  fact  that  Zq  =  u;/z0//?]0. 
Furthermore,  it  is  evident  that 


A'p  =  3(Zm) 


(269) 


-  u>/i06 


tan  kb\ 
kb  ) 


OO  OO 

E  E 

m  =  l  n=Q 
(”».*»)/  (1,0) 


en  sin2(m7rxo/a)e  °m"r 
amna  [1  -  ( wr/ kb )2] 


where 

<*mn  =  (mx/a)2  +  (nn/b)2  -  k2  (270) 

2 Recall  that  for  the  purpose  of  evaluating  Ztn  we  have  taken  the  guide  to  be  matched  at 
z  —  ±oo. 
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This  completes  the  formal  evaluation  of  the  parameters  necessary  to  specify  the 
equivalent  circuit  of  Figure  105.  We  now  reduce  the  expression  in  (269)  for  Xp  to 
a  simpler  and  more  easily  computed  form.  First  write  Xp  as 


X„ 


e  ainT  sin2(rxo/a) 
olna[l  —  ( nr /kb )2] 


(271) 


+  V - - - V 

i  1  -  (nr/ kb )2  “0 


tn  v" »  .  2  nirx0 

, — . .  >  sin  - 

( W M)2  a 


e-Omnf  g— nurr/o  ] 


rnr 


h 3  1  -  (mr/kb)2  ^  a  mrr  j 


In  what  follows,  we  set  x0  =  a/2;  this  is  the  situation  of  practical  interest.  Now 

[151 

y «  _ fn_ 

1  -  (nr, 


=  fc&cot 


„=o  -  -  (nr /kb)2 

y  J_e- nurr/a  =  _L  in  coth  -  IC~W« 
,  m7r  2r  \  2a)  r 


(272) 

(273) 


m=3 

(odd  m 


so  that 


(t&n  kb\  [  1  /  7rr\  e  ,rr/o 

Ap  =  uiLiob  — — —  )  <  — —  In  (  coth  — - 

p  A  \  Jfc6  J  \2r  V  2a J  7r 


(274) 


/tan  A:6\  ^ 

+  2(  -rr-l  E 


a,na  [1  -  (nr /kb)2] 


+ 


( tan  kb\ 

oo 

V  - 

Cn  y' 

g  —  Ofmn7" 

g  —  mnr/a' 

' 

V  kb  j 

^  1 
n=0  1 

-(nr /kb)2  ~ 

( odd  m) 

mr 

4 

The  infinite  series  in  this  representation  for  Xp  are  rapidly  convergent  and  the 
numerical  evaluation  of  Xp  is  straightforward.  Curves  of  A/,,  normalized  by  the 
free-space  intrinsic  impedance  r/0  =  \20r  fl,  are  shown  in  Figure  106  as  a  function 
of  normalized  frequency  ka/r  for  b/a  =  0.5  and  rja  =  0.01,  0.03,  and  0.1.  We 
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note  that  Xp  is  negative,  indicating  that  the  probe  impedance  is  capacitive  over 
this  frequency  range. 


ka/n 

Figure  106.  Normalized  probe  reactance  vs.  normalized  frequency  for  b/a  = 
1/2;  r/a  =  0.01,  0.03,  0.1. 


6.3  NUMERICAL  RESULTS  FOR  ABSORBED  POWER. 


In  this  section  we  present  numerical  results  for  the  normalized  absorbed  power 
of  Eq.  (253)  as  a  function  of  normalized  frequency  ka/n.  The  waveguide  and 
termination  dimensions  are  chosen  such  that 


b  =  a/2  (275) 

l  =  a/2  (276) 

This  choice  for  l  makes  Xj  =  oo  at  ka/n  =  \/2.  The  normalized  absorbed  power 
can  be  written 


Pats  =  _ \N2R<jZo _ 

Pine  (Rd  cot  0U>1  -  Xpy  +  (Rd  +  A2Z0  +  Xp  cot  AoO2 


(277) 
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where  R d  is  the  small-signal  diode  resistance,  Xp  is  given  in  (274),  N  is  given  in 
(268),  and 


(278) 

(279) 


We  have  performed  numerical  calculations  of  Paba/Pinc  as  a  function  of  kafir  for 
a  variety  of  cases.  First,  in  Figure  107  we  plot  Paba/Pinc  vs.  kafir  when  the  probe 
reactance  Xv  =  0.  In  this  case,  we  assume  that  the  reactance  is  tuned  out  at 
every  frequency;  Pabs/Pinc  becomes 


Figure  107. 


Pabjp, 


xp=o 


4N2RdZ0 

R%  cot2  (li0l  +  {Rd  +  N2Z0)2  ' 


1.0  1.2  1.4  1.6  1.8  2.0 

ka/n 


(280) 


Normalized  absorbed  power  vs.  normalized  frequency  for  Xp  =  0; 
Rdho  =  0.25,  0.495,  1.0. 
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This  quantity  has  a  maximum  with  respect  to  Rd  occurring  at 


Rd  =  N2  Z0  sm  PiqI 


(281) 


The  value  of  the  maximum  is 


_  2  sin  ftu,/ 

max;  1  -j-  sin  diQ l 

Xp=0 

When  /3i0l  =  7r/2,  a  condition  that  would  denote  the  design  point,  the  maximum 
value  of  Pabs/Pinc  is  unity.  A  value  of  Rd  =  0.495^0  meets  the  condition  of  Eq. 
(281)  at  the  design  point. 

We  observe  from  Figure  107  that  changing  Rd  from  this  value  still  permits  Pabs/Pinc 
to  approach  unity;  the  frequency  at  which  this  occurs  decreases  if  Rd  >0.495t7o 
and  increases  if  Rd  <  0.495r/o.  This  is  a  consequence  of  the  frequency  dependence 
of  N2Z0  and  of  cot  /310/.  We  conclude  from  this  result  that  if  the  probe  reactance 
can  be  tuned  out  (as  it  certainly  can,  at  least  at  a  given  frequency),  essentially  all 
of  the  power  incident  upon  the  diode  can  be  absorbed  by  it,  at  some  frequency  in 
the  band. 

We  show  curves  of  Pabs/Pinc  vs.  ka/ir  for  various  uncompensated  probe  induc¬ 
tances  (determined  by  differ  at  ratios  r/a)  and  various  values  of  Rd/r/o  in  Figures 
108  -  110.  The  trends  are  adequately  shown  by  the  curves  in  Figure  109  for  r/a 
=  0.03.  Pabs/Pinc  reaches  a  higher  maximum  value  at  a  lower  frequency  as  Rdho 
is  increased.  The  differences  in  peak  values  become  more  pronounced  as  r/a  de¬ 
creases,  that  is,  as  the  probe  inductance  increases.  However,  the  decrease  in  peak 
value  of  Pabs/Pinc  when  Rd/^o= 0.25  and  r/a  =  0.01  is  only  a  factor  of  two,  or 
3  dB.  We  conclude  that  even  when  the  probe  reactance  is  uncompensated,  the 
fractional  power  absorbed  in  the  load  can  be  substantial.  We  remark  also  that 
the  half-power  width  of  these  curves  is  nearly  ka/2ir  -  half  of  the  TEi0  operating 
band. 
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ka/n 

Figure  108.  Normalized  absorbed  power  vs.  normalized  frequency  for  r/a 
0.01;  Rj/tiq  =  0.25,  0.495,  1.0. 


Figure  109.  Normalized  absorbed  power  vs.  normalized  frequency  for  r/a 
0.03;  Ri/r) o  =  0.25,  0.495,  1.0. 
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ka/it 


Figure  110.  Normalized  absorbed  power  vs.  normalized  frequency  for  r/a  = 

0.1;  Rd/ %  =  0.25,  0.495,  1.0. 

6.4  CONCLUDING  REMARKS. 

We  have  constructed  a  mathematical  model  for  the  power  absorbed  in  a  resistive 
load  coupled  to  a  waveguide  by  a  thin  probe  extending  across  the  short  dimension 
of  the  guide  cross-section.  It  has  been  assumed  that  at  TEj0  mode  field  is  incident 
upon  the  probe  and  that  the  guide  is  terminated  in  a  short  circuit.  The  frequency 
range  considered  is  that  over  which  only  the  dominant  TEjo  mode  can  propagate. 
Representative  calculations  have  shown  that  even  when  the  probe  reactance  is 
uncompensated,  a  substantial  fraction  of  the  incident  power  is  absorbed  in  the 
load.  The  width  of  the  absorption  curve  is  approximately  one-half  of  the  TEjo 
mode  band  for  a  waveguide  whose  cross-sectional  aspect  ratio  is  2,  a  typical  value. 

The  values  of  Rd,  the  load  resistance,  which  were  considered  were  in  the  range 
100  -  400  fl.  This  range  includes  typical  small-signal  diode  resistances. 

It  is  not  surprising  that  coupling  to  a  load  in  a  waveguide  can  be  so  highly  efficient 
-  after  all,  this  is  the  desired  purpose  of  the  configuration.  From  the  standpoint 
of  vulnerability  modeling,  however,  the  relatively  broad-band  character  of  the 
response  over  the  frequency  range  associated  with  pure  TEjo  mode  propagation 
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would  lead  us  to  conclude  that  in-band  vulnerability  of  a  microwave  system  is 
largely  defined  by  the  power  which  reaches  the  diode  end  of  the  system.  Power 
which  reaches  this  point  will  be  absorbed. 

Time  has  not  permitted  us  to  explore  the  out-of-band  coupling  problem.  When 
more  than  one  mode  can  propagate  in  the  waveguide,  the  model  of  Figure  105  will 
contain  a  transformer  for  each  propagating  mode.  This  case  can  be  explored  using 
the  methods  employed  herein,  and  is  worthy  of  consideration.  Other  cases  which 
might  profitably  be  examined  include  that  in  which  the  probe  extends  only  part¬ 
way  into  the  waveguide  and  that  in  which  the  resistive  load  is  placed  inside  the 
waveguide  itself.  The  analysis  of  both  of  these  latter  problems  is  straightforward. 
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SECTION  7 

ELECTROMAGNETIC  PENETRATION  OF  LONG  SLOTS: 
AN  APPROXIMATE  ANALYSIS 


Summary:  Electromagnetic  penetration  of  an  infinitely  long,  electrically  narrow 
slot  of  arbitrary  depth  in  a  plane  conductor  is  investigated  using  a  relatively  simple 
analytical  technique.  The  quantity  of  interest  is  the  total  power  transmitted 
through  the  slot  under  plane-wave  illumination  with  magnetic  field  parallel  to 
the  long  axis  of  the  slot.  Analytical  and  numerical  results  for  this  quantity  are 
provided,  as  is  an  elementary  equivalent  circuit  model  of  the  problem.  The  analysis 
approach  can  be  applied  to  more  complicated  problems. 

There  exists  considerable  need  for  models  of  various  electromagnetic  coupling 
configurations  for  use  in  the  study  of  electromagnetic  interactions  with  systems  of 
many  kinds.  Unfortunately,  only  a  few  such  configurations  have  been  studied  in 
detail,  yet  the  need  for  models  remains.  At  the  present  time,  extensive  numerical 
and  experimental  programs  are  under  way  whose  purpose  is  to  develop  increased 
understanding  of  coupling  through  thick  apertures,  multiple  apertures,  loaded 
apertures,  penetrating  conductors,  and  other  complex  configurations.  While  these 
efforts  have  yielded  useful  results,  these  results  are  often  not  in  a  form  which  is 
readily  usable  by  one  responsible  for  making  approximate  coupling  estimates  or 
for  quantifying  hardening  requirements.  Analytical  approaches,  where  they  can 
be  applied,  are  often  the  best  means  for  obtaining  relatively  simple,  although 
sometimes  approximate,  results  for  near-term  use. 

Our  principal  interest  is  in  developing  relatively  simple  analytical  expressions  for 
the  power  transmitted  through  a  rectangular  slot  of  arbitrary  depth,  illuminated 
by  a  plane  wave.  In  this  section  we  discuss  the  basis  for  our  approach  to  that 
problem  in  the  context  of  a  simpler  one:  the  case  in  which  the  slot  is  infinitely 
long.  This  problem  has  already  been  treated  by  Harrington  and  Auckland  [16] 
using  a  slightly  different  method  from  that  employed  herein.  The  analysis  and 
the  results  developed  in  this  chapter  are  not  new;  but  it  is  our  intent  to  illustrate 
this  approach  for  treating  other  more  complicated  problems. 

The  method  which  we  shall  employ  herein  is  described  in  Harrington’s  text  Time- 
Harmonic  Electromagnetic  Fields  [17].  It  consists  in  making  reasonable  assump¬ 
tions  concerning  the  tangential  electric  field  in  the  aperture  connecting  two  regions 
of  space,  calculating  the  electromagnetic  field  in  each  region  from  the  assumed 


141 


aperture  field,  and  then  requiring  that  the  complex  power  flow  be  continuous 
through  the  aperture.  This  final  condition,  which  permits  the  solution  for  the 
aperture-field  amplitude  to  be  obtained,  represents  an  approximate  means  for 
satisfying  the  condition  that  the  tangential  magnetic  field  be  continuous  through 
the  aperture.  This  technique  is  ideally  suited  to  many  of  the  problems  encountered 
in  microwave  coupling  to  systems  via  back-door  paths. 

We  begin  by  considering  the  problem  of  penetration  of  a  long,  narrow  slot  in  an 
infinitesimally  thin  conductor.  This  serves  to  illustrate  the  method,  to  demon¬ 
strate  that  the  results  are  not  sensitive  to  the  form  of  aperture  field  chosen,  and 
to  show  how  a  simple  equivalent  circuit  can  be  derived  which  contains  all  the 
relevant  features  of  the  original  problem.  Then  we  address  the  deep-slot  problem. 
We  construct  an  equivalent  circuit  and  present  representative  numerical  results 
indicating,  inter  alia ,  the  slot-depth  resonance.  The  formulation  of  this  problem 
includes  the  possibility  of  slot  loading  by  a  material  different  from  free  space.  We 
also  present  a  brief  discussion  of  the  problem  of  coupling  to  an  infinitely  deep  slot. 
Certain  mathematical  results  are  given  in  the  Appendix. 

7.1  THE  THIN,  LONG  SLOT:  FORMULATION  AND  SOLUTION. 

The  geometry  of  the  problem  is  shown  in  Figure  111.  An  H-polarized  plane 
electromagnetic  wave  is  incident  at  an  angle  0  (|  0  |<  tt/2)  on  a  thin,  perfectly 
conducting  plane  occupying  the  surface  2  =  0.  The  plane  is  cut  by  an  infinitely 
long  slot  aperture  in  the  region  |  x  |<  iu/2;  the  center  of  the  slot  lies  along  the 
j/-axis.  The  medium  on  either  side  of  the  conducting  plane  is  free  space.  It  is 
assumed  that  the  slot  is  electrically  thin;  that  is,  for  field  time  variations  of  the 
form  exp  {jut),  kw  =  unv/c  <<  1  where  c  denotes  the  speed  of  light.  The  object 
of  the  analysis  is  to  calculate  the  power  transmitted  through  the  slot. 

Our  approximate  approach  to  the  solution  of  this  boundary- value  problem  pro¬ 
ceeds  as  follows:  first,  an  electric  field  Ex(x)  of  (assumed)  known  spatial  variation 
but  unknown  amplitude  is  imposed  across  the  slot  and  the  total  electromagnetic 
field  is  computed  in  the  regions  z  <  0  and  z  >  0.  Next,  the  fields  are  connected 
through  the  slot  by  forcing  the  complex  power  flow  to  be  continuous  there.  This 
condition  permits  us  to  obtain  the  amplitude  of  the  electric  field  across  the  slot 
and  thus  to  calculate  the  power  transmitted  through  the  ’ot. 


142 


The  electromagnetic  field  comprises  components  Ex,  Hy,  and  Ez  related  by 

Ex  = 


where 


E,  = 


\dx2 


(u o *co> 


-ZodHy 
jk  dz 

(283) 

ZodHz 
jk  dx 

(284) 

■  +  jfc2  j  =  0 

(285) 

(Vco) 

W 

Figure  111.  Geometry  of  the  thin-slot  problem.  The  slot  width  w  is  assumed  to 
be  small  compared  to  the  wavelength. 
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and  Zq  denotes  the  intrinsic  impedance  of  free  space.  Denoting  by  E0  the  ampli¬ 
tude  of  the  incident  electric  field  and  by  V0  the  total  voltage  across  the  slot 


Vo  =  - 


(286) 


we  obtain  the  following  expressions  for  the  relevant  components  of  the  electro¬ 
magnetic  field: 


z  <  0  :  Hy 


~e-jkxsin6cos(kzcosd) 

Zq 


+ 


kVQ 

2irZ0 


jhx+js/k1  —h2z 


dh 

Vk2  -  h2 


(287) 


Ex  =  2jE0co$Oe  2kxsme  sm(kz  cosO) 


(288) 


K> 

27T 


j  hx+j  '/k2-h2  z 


dh 


z>0:Hy~  2 xZ 


-kv0  r 


Ex  = 


r  r  F(h)e}hx-J'/kJ^z—/  d  = 

o  J-oo  \/k2  —  h2 

r  F(h)<zjhx->Vki-hizdh 

J  —  OO 


-Vo  f°° 

27T 


(289) 

(290) 


The  function  F(h)  is  to  be  determined.  The  field  expressions  given  above  are  read¬ 
ily  shown  to  satisfy  Eqs.  (283)  and  (285)  as  well  as  proper  boundary  conditions 
at  |  2  |— >  oo  if  _ 


KVV-h2  >  0 
%Vk2  -  h2  <  0 


(291) 

(292) 


The  electric  field  in  the  aperture  plane  Ex(x,0)  is 

Ex(x,0)  =  ^  r  E(li)ejhxdh 

2ir  J—oo 


whence 

F(h)  =  —77  rn  E,{i',0)e-ihx‘dx‘ 
Vq  J—w/2 


(293) 

(294) 
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and  from  Eq.  (286),  F( 0)  =  1.  Thus  the  tangential  electric  field  Ex  is  continuous 
at  z  =  0  and  Ex  =  0  on  the  conducting  surface  if  Eq.  (294)  is  satisfied. 

An  exact  solution  to  this  problem  would  follow  from  forcing  continuity  of  Hy 
through  the  slot.  This  condition  can  be  written 


—ikx  sin  9 

Z0 


-kV0 

■kZq 


dh 

s/k 2  -  h 2 


(|  X  |<  w/2) 


(295) 


which  can  easily  be  shown  to  lead  to  an  integral  equation  for  the  electric  field 
across  the  slot.  We  shall  apply  a  less  stringent  condition:  multiply  both  sides 
of  Eq.  (295)  by  (l/2)i?*(:r,0)  and  integrate  over  all  x.  This  condition  makes 
the  complex  power  continuous  through  the  slot,  a  necessary  but  not  sufficient 
condition  for  the  satisfaction  of  Eq.  (295).  We  obtain 


1 

2 


0-jkx  sin  6 


5  £  *  (ir )  £  rw*-,k’dh 


(296) 


dh' 


v/F  -  h '2 


and  making  use  of  the  result 


e-jo<d(,  =  6(a) 


(297) 


where  <5(-)  denotes  the  Dirac  delta-function,  we  find  that  Eq.  (296)  reduces  to 

~kVo  r°  ,  12  dh 


E0F*(—k sin 0)  =  /°°  |  F(/l) 

7T  JO 


n/F  -  h 2 

Furthermore,  the  real  power  transmitted  through  the  slot  is  given  by 


R  = 


k  |  Vo  |2  [k 


I  \m 

Jo 


dh 


v/F  -  h2 


2ttZo 

Now  it  is  evident  from  Eq.  (294)  that 

F(  —  ksin0)  =  1  +  £?(&tysin  0)2 

Further,  defining 


roo 

/  I  F(h) 
Jo 


dh 


Vk2  -  h 2 


—  It  +  jU 


(298) 

(299) 

(300) 

(301) 
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we  have 


7'  =  i|FW|!T TTT^l  (tw<<1> 

Ii  =  J°°  I  /•(*)  I2  rfA 


so  that,  for  kw  <<  1, 


and 


K>  = 

Pt  = 


Vh*  -  fc2 
—2Eqw 

kw(l  +  2jli/n) 
w\  E0\2 


Z0kw  [l  +  (2/,/tt)2] 
We  may  define  an  effective  slot  width  we  as 

^inr 


(302) 

(303) 

(304) 

(305) 

(306) 


where  5;nc  =|  E0  |2  /(2 Z0)  denotes  the  power  density  incident  on  the  slotted  plane. 
We  obtain 

2  w 

We  =  kw  [1  -f  (2/,-/ 7T )2]  (307) 

whence 

We  _  _ 2 _ 

w  kw  [1  +  (2/,/7t)2] 
is  the  normalized  effective  slot  width. 


(308) 


It  is  useful  to  construct  an  equivalent  circuit  describing  the  penetration  of  the 
slot.  From  Eq.  (298)  with  F*(  —  fcsinfl)  =  1,  we  have 

Vo  =  -J;;/(2Y')  (309) 

where  Jascx  =  2£0/Z0  is  the  short-circuit  current  density  across  the  slot  and 

k 


is  an  admittance  per  unit  length.  Writing  Y'  =  G'  +  jB',  we  have 


G'  =  ^z 


(310) 


(311) 
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B'-bu  (312) 

where  G'  and  B'  are  respectively  the  slot  conductance  and  susceptance  per  unit 
length.  The  equivalent  circuit  is  shown  in  Figure  112. 

It  remains  to  consider  the  slot  electric  field  and  evaluate  the  integral  To  show 
that  the  result  is  not  sensitive  to  the  specific  form  chosen  for  Ex(x,  0)  we  consider 
two  different  expressions  for  this  field: 


£<1}(*,0) 


(313) 


El2)(x,0) 


(314) 


The  first  of  these  expressions  represents  a  field  which  is  constant  across  the  slot. 
From  Eq.  (294),  we  obtain 


F(h)  ~  sine 


(315) 


where  sinc(x)  =  sinx/x.  The  second  expression  for  Ex(x,0)  accounts  for  the 
singular  behavior  of  the  electric  field  at  the  edges  of  the  slot  and  is  an  accurate 
representation  of  the  field  in  the  limit  kw  — >  0  [18].  Again  from  Eq.  (294),  for 
this  case 

F(h)  =  J„  (316) 

where  J0(-)  is  the  Bessel  function  of  order  zero.  Thus  define 


-  f, 


00  sine2  (kwu/ 2)  du 
\Ju 2  —  1 

j( 2)  _  yT°°  Jp  (kiuu/2)  du 
It  is  shown  in  the  Appendix  that 


y/u2  —  \ 


lim  Z,'1’  =  3/2  +  In  (  — — 

kw^O  l  -y  ktV 


(317) 

(318) 


(319) 
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where  7  =  1.781-  ■  •  is  the  exponential  of  Euler’s  constant.  It  is  apparent  that  the 
change  produced  in  7,  by  different  choices  for  the  slot  electric  field  occurs  only  as 
a  factor  multiplying  lew  in  the  argument  of  a  natural  logarithm,  a  slowly  varying 
function. 

We  show  in  Figure  113  plots  of  ive/iu  as  a  function  of  kiv/r  for  7,  =  7-1'  and 
Ii  =  7,-2*.  One  will  observe  that  the  differences  are  small,  a  result  of  the  fact  that 
the  quantity  calculated  depends  on  the  integral  of  the  aperture  field  rather  than 
on  its  detailed  behavior.  We  shall  use  this  fact  to  advantage  in  the  next  section, 
wherein  we  consider  the  effect  of  finite  slot  depth. 


(illuminated  side) 


(shadowed  side) 


Figure  112.  Equivalent  circuit  for  the  thin-slot  problem.  J is  the  short-circuit 
current  density  across  the  slot  and  Y'  is  an  admittance  per  unit 
length. 


Figure  113.  Normalized  effective  slot  width  we/w  as  a  function  of  kiv/x. 

The  solid  curve  results  from  choosing  U  —  I,-1*  and  the  dashed 
curve  from  /,  =  lj2\  xf 


I— d— I 


Figure  114.  Geometry  of  the  deep-slot  problem.  It  is  assumed  that  kw  <<  1 
and  that  kd  is  arbitrary.  The  slot  may  be  filled  with  a 
homogeneous,  isotropic  medium  other  than  free  space. 


7.2  THE  DEEP,  LONG  SLOT:  FORMULATION  AND  SOLUTION. 


The  geometry  of  the  deep-slot  problem  is  shown  in  Figure  114.  The  depth  of  the 
slot  is  d  and  it  is  assumed  that  this  depth  is  not  necessarily  small.  As  in  the  thin- 
slot  problem,  we  shall  assume  that  the  forms,  but  not  the  complex  amplitudes,  of 
the  electric  fields  across  the  entrance  and  exit  of  the  slot  are  known. 


The  fields  in  the  regions  outside  the  slot  are  given  by 
2  <  0: 


H  -  2Eo 

Hy  - 

^0 


-jkx  sin  6 


cos (kz  cos  0) 


+ 


kVj 
2t rZ, 


/OO 

F(h), 

■OO 


Jkx+jVk2  -/»2  z 


dh 


s/l t2  -  h 2 


-2jE0  cos  0e-}kxsine  s\n(fcz  cos  6) 


(321) 


(322) 


z  >  d: 


_K 

2tt 


j/ix+ j \/ k2~~h?  z 


dh 


Hy  2  tZ, 


/  /■OO 

~  /  F(/i)ejfcr~J 

0  ^ — oo 


tx—ji  \/fc*  —h?  { z—  d )  _ 


x/A:2  -  h 2 


Ex  =  — ^  [°°  F(h)e3hz-lVkTZJ?{!-d)dh 

2n  7-oo 


(323) 

(324) 


where  now  V;  denotes  the  voltage  across  the  input  aperture  at  2  =  0  and  Vo  is  the 
voltage  across  the  output  aperture  at  z  =  d.  We  shall  assume  that  the  slot  field 
is  uniform,  so  that 


F(h)  =  sine 


(325) 


The  field  in  the  slot  itself  is  represented  by  forward-and  backward-going  TEM 
modes.  In  order  to  account  for  loading  within  the  slot,  we  assume  that  the  prop¬ 
agation  constant  k,  and  intrinsic  impedance  Zs  of  the  medium  filling  the  slot  are 
not  necessarily  those  of  free  space.  Thus 
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0  <  z  <  d: 


Ex  =  Ae~jk’z  +  Bejk’z  (|  x  |<  u>/2) 

Hy  =  ~  (Ae~jk‘z  -  Be]k‘z)  (|  x  |<  w/2) 

where  the  amplitudes  A  and  B  are  determined  from  the  relations 

A+  B=~ 
w 


and  are  given  by 


Ae~jk’d  +  Be]k’d  =  — - 
w 

A  =  (2jiw  sin  k,d)-'(-V,ejk’d  +  V0) 
B  =  (2jie  sin  kad)~l(  —  Vo  +  Vi^k,d) 


(326) 

(327) 


(328) 

(329) 

(330) 

(331) 


As  in  the  thin-slot  problem  considered  in  the  previous  section,  we  construct  equa¬ 
tions  expressing  the  continuity  of  Hy  across  the  slot  apertures  and  then  satisfy 
these  conditions  in  an  integral  sense.  Continuity  of  Hy  at  z  =  0  yields  the  condi¬ 
tion 


-jkx sing  ,  W 

Z0 


k  V  r°° 


jhx 


dh 


Vk2  -  h 2 


(332) 


1 


jwZa  sin  ksd 

Continuity  of  Hy  at  z  —  d  yields 


(-Vi  cos  kad  +  V0)  ((  x  |<  u>/2) 


- kVo 

2k  Z{ 


'c\  r°° 

r  /  F(h)elhx- 
0  J-OO 


dh 


1 


\A2  —  h2  ju>Za  sin  kad 


(-Vi  +  V0  cos  kad)  (|  x  i<  tu/2) 

(333) 

Now  multiply  Eq.  (332)  through  by  (1/2)  E*(x,0)  and  Eq.  (333)  by  (1/2 )E^.(x,d) 
and  integrate  with  respect  to  x.  The  result  is  the  pair  of  equations 


*  +  f  (■+¥) 

Zq 


(334) 


2 jwZa  sin  k,d 


(V0- Vi  cos  kad) 
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(335) 


(l  +  M),  .h  (VoeosM-  Vi)  (335) 

4  \  ir  )  2jwZa  sin  kad 

where  we  have  used  the  fact  that  for  kw  «  1,  F*(—ksmO)  =  1  and  Ir  =  ir/2. 

Now  it  is  a  simple  matter  to  solve  Eqs.  (334)  and  (335)  for  the  slot  voltages  V{ 
and  V0.  We  obtain 


E0w  [(Z0/2jZa)cot  kad  +  (kw/ 4)  (l  +  ^/.)] 


E0w  [( ZQj  2j  Za )  esc  ksd\ 


where 


esc  kad\  — 


Zo  .  ,  ,  ,  kxv 

W.cot  m+T 


(1  +  l/‘)]! 


(336) 


(337) 


(338) 


As  in  the  thin-slot  case,  we  can  construct  an  equivalent  circuit  representing  the 
thick-slot  configuration.  This  equivalent  circuit  is  shown  in  Figure  115.  The 
admittance  Y  is  given  by 

K  =  ^(1  +  T/')  <339> 


and  the  source  current  /,  is 


Is  =  r>  (340) 

This  equivalent  circuit  yields  the  correct  voltages  at  the  slot  apertures  and  at 
points  inside  the  slot.  Powers  calculated  with  this  circuit  should  be  divided  by 
the  slot  width  w  to  yield  power  per  unit  length  along  the  slot. 

The  normalized  effective  width  of  the  slot  is 


wjw  = 


kw\  Zq  esc  kad 
T  2  jZaD 


(341) 
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Figure  115.  Equivalent  circuit  for  the  deep-slot  problem. 


which  can  be  written  as 


we/w  = 


\j  (Z0/Z,  +  Y2Z0Z,)  sin  kad  +  2 V Z0  cos  kad\' 


(342) 


Note  that 


lim  w./w 
d~*  o 


kw  1  +  (2/,/7t)2 


(343) 


as  required  by  Eq.  (308). 


Curves  of  we(d)/we(0)  are  plotted  in  Figure  116  as  a  function  of  kd/ir  for  fixed 
values  of  kiv/ir  when  ka  =  k  and  Z,  =  Zo,  i.e.,  when  the  slot  is  unloaded.  We 
observe  an  initial  decrease  in  we(d)/we( 0)  as  kd/ir  is  increased  from  zero,  followed 
by  an  increcise  when  kd/ir  is  close  to,  but  less  than,  one,  and  return  to  unity 
at  kd/ir  =  1.  This  effect  is  readilv  understood  in  terms  of  the  circuit  model  of 
Figure  115.  As  kd/ir  increases,  the  admittance  seen  looking  into  the  entrance 
aperture  changes.  Between  kd/ir  =  0  and  1  there  will  exist  a  condition  where  a 
conjugate  match  of  source  and  apparent  load  admittances  will  most  nearly  occur; 
the  transmitted  power  will  be  maximized  at  this  point.  The  maximum  possible 
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values  of  we/w  and  we(d)/we( 0)  can  be  found  from  the  maximum  power-transfer 
theorem  to  be 


we/w 
we{d)/we(  0) 


2 

kw 

(344) 

II 

+ 

W 

(345) 

It  is  of  potential  practical  interest  in  revealing  the  existence  of  deep-slot  penetra¬ 
tions  to  consider  the  effective  width  of  the  entrance  aperture  for  power  scattered 
back  into  the  illuminated  region.  This  normalized  effective  width,  denoted  w'e/w , 
is  given  by 


w'e/w  =  ^ 

Curves  of  w'e(d)/we{ 0)  are  plotted  as  a  function  of  kd/ir  for  various  values  of  lew/ n 
and  k ,  =  k ,  Z,  =  ZQ  in  Figure  117.  We  observe  that  w'e (d)/we(0)  behaves  similarly 
to  we(d)/w( 0).  Thus  it  would  appear  that  deep-slot  resonance  effects  should  be 
observable  from  either  side  of  the  slot. 

7.3  PENETRATION  OF  AN  INFINITELY  DEEP  SLOT. 


kw' 

T, 


(Z0/2jZs)  cot  kad  +  (kw/ 4)  [1  +  (2jlj/n)\ 
D 


(346) 


If  the  slot  is  taken  to  be  infinitely  deep,  the  fields  in  the  slot  are  of  the  form 


Ex  =  —  C]k-Z  (|  x  |<  w/2) 
w 

Continuity  of  Hy  across  the  opening  of  the  slot  requires  that 


(347) 

(348) 


-)fcr  sing  , 

Z0  2irZ( 


/oo 

F{h)elhx • 

■OO  1 


dh _ -Vj 

y/k2  -  h2  ~  Z,w 


(|  x  |<  w/2)  (349) 
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Figure  116.  Normalized  effective  width  we(d) /ive(Q)  as  a  function  of  kd/x  for 
kw/ir  =  0.01,  0.03 ,  0.1.  The  slot  is  unloaded,  i.e.,  k ,  =  k 
and  Z3  —  Z0. 


Figure  117.  Normalized  effective  width  w'(d)/we( 0)  as  a  function  of 

kd/x  for  kw/w  =  0.01,  0.03,  0.1.  The  slot  is  unloaded,  i.e., 
kf  =  k  and  Zs  —  Zq. 
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Enforcing  this  condition  in  the  integral  sense  used  throughout  this  section  yields 
an  expression  for  K: 


—2  wEq 

Zo/Z,  +  (kw/2)(l+2jIi/7r) 


The  real  power  entering  the  slot  is 


and  the  normalized  effective  slot  width  is 

4Z0dl(l/Za) 

W,W  |  Zo/Z,  +  ( W  2)  (1  +  2j/,/ir)|2 


(350) 


(351) 


(352) 


A  plot  of  we/w  as  a  function  of  kw/ir  is  shown  in  Figure  118  for  the  case  Za  =  Z0, 
i.e.,  when  the  slot  is  empty.  Note  that 


lim  we/w  =  4%t(Za)/Z0 

kw— >0 


(353) 


Figure  118.  Normalized  effective  area  of  an  infinitely  deep  slot  as  a  function  of 
kw/n.  The  slot  is  taken  to  be  unloaded,  i.e.,  ka  =  k  and 
Za  =  Zq. 
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7.4  CONCLUDING  REMARKS. 


We  have  presented  some  approximate  analysis,  equivalent  circuits,  and  relatively 
simple  formulas  for  assessing  the  penetration  of  electromagnetic  fields  through 
slot  apertures  which  are  electrically  narrow  (the  short  aperture  dimension  is  small 
compared  to  the  wavelength),  infinitely  long,  and  of  arbitrary  depth.  We  have 
employed  a  technique  which  can  be  utilized  in  many  other  related  problems  of  in¬ 
terest  in  microwave  coupling,  including  penetration  through  a  slot  of  finite  length. 

We  wish  to  stress  in  concluding  this  chapter  that,  while  not  exact,  analytical  meth¬ 
ods  such  as  those  employed  herein  can  provide  useful  models  for  many  “back-door” 
coupling  problems.  They  should  be  “bench-marked”  by  checking  the  limiting  cases 
for  which  exact  analytical  solutions  are  available  and/or  by  careful  numerical  mod¬ 
eling  studies  of  special  cases.  Through  this  means  the  limitations  on  the  validity 
of  the  simple  analytical  models  can  be  established  and  their  accuracy  quantified. 
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APPENDIX  TO  SECTION  7 
EVALUATION  OF  INTEGRALS 


EVALUATION  OF  ljl)  AS  kw  -+  0 


We  express  sine  (kivu/2)  in  integral  form  as 


sine 


(kw  \ 

1 

1  — tl 

—  ~ 

V  2  i 

2  7-, 

Jkwu, 


(354) 


and  write  /-°  as 


/'"  =  i  /'  f  r  (355) 

4  7-i7-i  7i  y/u2  —  1 

=  T  /,  (t  ' f  ” "  •) 

where  Vo(-)  denotes  the  Neumann  function  of  order  zero.  Since  kw  <<  1,  the 
argument  of  the  Neumann  function  will  be  small  over  the  range  of  integration  and 
the  function  can  be  replaced  by  its  small-argument  approximation 


lim  Y0{x)  =  —  In  — 


r— *0 


,7X . 


(356) 


where  7  =  1.7S1  ■  •  •  is  the  exponential  of  Euler’s  constant.  Thus 


lim  /,(,)  =  -  f  I'  hi  f  — - - - -]  d^h) 

kw—o  4  7-1  7_i  \7Aur  j  C  —  77 


which  is  readily  evaluated  to  yield 


...  3  9 

lim/,'0  =  -  +  ln— - 

k~o  '  2  7 kw 


(357) 


(358) 


=  In 


8.963’ 
7  k  w 
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EVALUATION  OF  lj2)  as  lew  ->  0 

We  use  the  integral  representation  of  the  Bessel  function 


Jo  (t“)  =  h  f_/kwucos4,/2d<t> 


(359) 


to  write  /-2*F  in  the  form 


I 


(2) 


=  jl  r  r  dm  r ~sikwu  060) 

(27T)  J~  -r  J—n  Jl  u?  —  1 

=  /  Y°  (~v  I  cos ^  -  cos 0  I  j  d<t>d<p 


Again  using  the  small-argument  approximation  for  the  Neumann  function,  we 
obtain 


lim  ll2)  =  f  r  In  f  — 

kw—0  47T2  J-  r  J-7T  \7«- 

and  with  the  help  of  the  identity 


w  |  cos  <j>  —  cos  i/>  |  / 


dcj>d^  (361) 


.  I 

In  (2  |  cos  <f>  —  cos  x!>  | )  =  —  2  }  —  cos  n<£  cos  nip 


n=l 


we  find 


lim  /,-2*  =  In  — - — 
fcu/— o 


(362) 


(363) 


In  - - (-  In  4 

7  Ktl) 


=  1.3863  + In 


7  A- in 
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